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Abstract. Wc consider a model case for a strictly convex domain 51 C M.'' of dimension 
d > 2 with smooth boundary dCl ^ and wc describe dispersion for the wave equation 
with Dirichlet boundary conditions. More specifically, we obtain the optimal fixed time 
decay rate for the smoothed out Green function: a t^^^ loss occurs with respect to the 
boundary less case, due to repeated occurrences of swallowtail type singularities in the 
wave front set. 

1. Introduction 

Let us consider solutions of the linear wave equation on a manifold {Q, g), with (possibly 
empty) boundary dfl: 

( {d^ - Ag)u{t, x) = 0, xen 
(1.1) < u{0,x) = uo{x), dtu{0,x) = ui{x), 

[ u(t, x) = 0, X e dfl , 

where denotes the Laplace-Beltrami operator on Q. 

When dealing with the Cauchy problem for nonlinear wave equations, one starts with 
perturbative techniques faces the difficulty of controlling the size of solutions to the linear 
equation in terms of the size of the initial data. Of course, one has to quantify this notion 
of size by specifying a suitable (space-time) norm. It turns out that, especially at low 
regularities, mixed norms of type L^L^ are particularly useful. Moreover, the arguments 
leading to such estimates turn out to be useful when considering spectral cluster estimates, 
which are of independent interest (see |15j). 

On any smooth Riemannian manifold without boundary, the following set of so-called 
Strichartz estimates holds for solutions of the wave equation (11.11) (for T < oo) 

(1-2) \\u\\Li{o,T)L^-{n) < Ct{ ||%|Ih/3(q) + ll^i|li//3-i(n)) , 

where, if d denotes the dimension of the manifold, we have /3 = d{^ ~ r) ^ g (which is 
consistent with scaling) and where the pair (g, r) is wave-admissible, i.e. 

(1.3) g>2, - + ^^—l<^!—l (g>2ifrf = 3andg>4if rf = 2). 
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When equality holds in (11. 3p we say that the pair (g, r) is sharp wave- admissible. Here 
denotes the Sobolev space over VL. If this holds for T = oo we say there are global 
Strichartz estimates. Such inequalities were long ago established for Minkowski space (fiat 
metrics) and can be generalized to any (fi, g) where g is smooth, at least locally in time 
(thanks to the finite speed of propagation). 

The canonical path leading to such Strichartz estimates is to obtain a stronger, fixed 
time, dispersion estimate, which is then combined with energy conservation, interpolation 

and TT* arguments to obtain (II. 2p . Let us denote by e^**v^~^^ the half- wave propagators 
in fiat space, and G C^(]0, oo\). The following dispersion inequality holds: 

(1.4) ||^(_/,2Aj,,)e±^*v^||^,(^,)^^^(i,,) < C(rf)/i-^min{l, (Vl^l)"^}- 

Our aim in the present paper is to obtain these estimates inside domains. In fact, [9] 
outlines a roadmap to prove such a dispersion estimate, on a finite time interval, for 
solutions of (II. ip inside a strictly convex domain {VL,g) of dimension d>2. A complete 
description of the geometry of the (semi-classical) wave front set is provided for the solution 
to (II. ip with initial data (uq, Ui) = {6a, 0), where a e is a point sufficiently close to the 
boundary (depending on the scale h). This wave front set has caustics developing in 
arbitrarily small times and this induces a loss of 1/4 in (ll.4p for the h/\t\ factor. 

In the present work, we aim at completing the roadmap by constructing a suitable 
parametrix for such a solution and then proving dispersion for the approximated solution. 
It should be noted that parametrices have been available for the boundary value problem 
for a long time, see [HI [TOl H] , as a crucial tool to establish propagation of singularities 
for the wave equation on domains. However, while efficient at proving that singularities 
travel along the (generalized) bi-characteristic flow, they do not seem strong enough to 
obtain dispersion, at least in the presence of gliding rays. In the outside of a strictly con- 
vex obstacle (no gliding rays), the Melrose- Taylor parametrix was utilized in [H] to prove 
Strichartz estimates hold as in the case. All other positive results ([2] and references 
therein) rely instead on reflecting the metric across the boundary and considering a bound- 
ary less manifold with a Lipschitz metric across an interface, and then using the machinery 
originally developed for low regularity metrics [131 [H] and spectral cluster estimates [15] . 
Such constructions do away with multiply reflected rays by suitable microlocalizations: one 
ends up working on a possibly very small time interval, depending on the incidence of the 
wave packet under consideration, such that all corresponding rays are only reflected once. 
Summing these intervals induces (scale-invariant) losses, which get worse with dimension; 
while Strichartz estimates are obtained in a more direct way in [2], one can observe that 
the corresponding dispersion estimate would have at most 1/t decay for c? > 4, as the ar- 
gument is blind to the full dispersion which should occur in tangential directions. On the 
other hand, negative results were obtained in [6l[5], where a special solution is constructed, 
propagating a cusp across multiple reflections and providing a counterexample to the sharp 
Strichartz estimates (II. 2p . for r > 4. This special solution is constructed via a microlocal 
parametrix which utilizes the Melrose one, and our present construction generalizes this 
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Figure 1.1. Light rays in tlie Friedlander model 
special example while retaining most of its useful features. 

Before stating our main result, we briefly introduce the Friedlander's model domain of 
the half-space = {{x, y)\x > 0,y & M.'^"^} with Laplace operator given by 

A, = dl + il + x)Ay. 

By rotational symmetry, we will eventually reduce to the two-dimensional case 

Remark 1.1. For the metric g = dx"^ + (1 + x)^^dy'^, the Laplace-Beltrami operator is 
Agfi = (1 + xY^'^dx{l + x)~^^'^dx + (1 + x)Ay which is self adjoint with the volume form 
y^detg dxdy = {1 + x)~^^'^dxdy. The Friedlander's model uses instead the Laplace operator 
associated to the Dirichlet form J iVgul"^ dxdy = J {\dxu\'^ -|- (1 -|- x)\dyu\'^)dxdy and is self 
adjoint with volume form dxdy. As a model, the Friedlander operator is better than the 
Laplacian A^ o since it allows explicit computations. Clearly, manifold [Qa^g) is a strictly 
convex domain : In fact, on the geodesic flow starting at x = 0,y = yo, ^0+^70 — g]0, 1 [, 
one has x{s) = 2s^q — s'^tjq. Moreover, (0,2,9) be seen as a simplified model for the 
disk -D(0, 1) with polar coordinates (r, 6), where r = 1 — x/2 and 6 = y. Multiply reflected 
light rays become periodic curves in the y variable, as illustrated on Figure 11.11 

Remark 1.2. We will always work with the Dirichlet boundary condition. The Neumann 
boundary condition can be handled exactly in the same way, providing the same results: 
one simply modifies the reflexion coefficient in our parametrix construction and replaces 
zeros of the Airy function Ai by zeros of its derivative Ai'. 

We are now in a position to state our main result. 

Theorem 1.3. Let d > 2. There exists C > 0,To > such that for every a g]0, 1], 
h G (0, 1] and t G (0,To], the solution Ua to (11. ip with data {uq,Ui) = (5a, 0) where 6a is 
the Dirac mass at point (a, 0, ■ ■ ■ ,0) G Qd, satisfies 

(1.5) \^i-h^Ag)ua{t,x)\ < Ch-'^ mm{l, {h/t)'^+i). 

The dispersion estimate (11.51) may be compared to (II. 4p : we notice a 1/4 loss in the h/t 
exponent, which we may informally relate to the presence of caustics in arbitrarily small 
times if a is small. Such caustics occur because optical rays are no longer diverging from 
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each other in the normal direction, where less dispersion occurs as compared to the M"^ 
case. We will prove in fact a slightly better estimate than fll.Sp : the {h/tY^^ factor may 
be replace by h^/^ + {h/tf/^ for a < h^'^ ( proposition [33D and by {h/ty/^ + a^/^h^'^ for 
a > h^^'^~^ ( Theorem 12 .11) . In fact, we can track the caustics and therefore our estimate is 
optimal for a > h^/"^^^ . 

Theorem 1.4. Let d > 2 and Ua he the solution to f 1 1.1 1) with data (mo,Mi) = (^ajO). Let 

h E (0, 1] and a > h^^'^~^ . There exists a finite sequence {tn)n, ^ ^ n < min(a^"'^/^, a^^'^h"^^^) 
such that 

(1.6) h-'^{h/tn)'^n-^/'^ah^/'^ < \ip{-h'^Ag)uaitr„a)\ and t^ ~ 4nv^, 

As a byproduct, we get that even for t g]0, Tq] with Tq small, the 1/4 loss is unavoidable 
for a small enough. We will see soon that this optimal loss is due to swallowtail type 
singularities in the wave front set of Ua- 

Remark 1.5. It follows from our proof that theorem 11.31 and 11.41 holds true if one replace 

^lj{-h^Ag)ua{t,x) by tlj{hDt)e^''V^^6,=a,y=o with V e C^{R*). 

As a consequence of (11. 5p and classical arguments, we obtain the following set of Strichartz 
estimates. 

Theorem 1.6. Let u be a solution of (II. ip on the model domain Qd, d > 2. Then there 
exists T such that 

(1-7) \\u\\Li{0,T)L^{n) < Ct{ \ \Uo\\hi3{Q) + ll'"lllH/9-i(n)) ) 

for {d, q, r) satisfying 

1 ,d-2 11 1 

and P is dictated by scaling. 

In dimension d = 2 the known range of admissible indices for which sharp Strichartz 
hold is in fact slightly larger, see [2]. However, in larger dimensions d > 3, Theorem 11.31 
improves the range of indices for which sharp Strichartz do hold, and does so in a uniform 
way with respect to dimension, in contrast to [2|. On the other hand, our results are, for 
now, restricted to a model case of strictly convex domain, while [2] applies to any domain. 
One may use the model case analysis to extend estimates to any smooth strictly convex 
domain, as in the counterexample situation [3]. This issue will be addressed elsewhere. 

Remark 1.7. One conjectures that the loss in Strichartz estimates in [6] are optimal. This 
would heuristically match a 1/6 loss in the dispersion estimate. We plan to address this 
issue in future work, by proving that the worst time-space points (t„, a) may be suitably 
averaged over. 

One may then make good use of such Strichartz estimates for the local (and global) 
Cauchy theory of nonlinear wave equations. We provide one simple example. 

4 

Theorem 1.8. The energy critical wave equation OgU + \u\d^u = with data {uq,Ui) G 
Hq^Q^) X L'^i^d) has unique global in time solutions for 3 < d < 6. 



DISPERSION FOR THE WAVE EQUATION INSIDE STRICTLY CONVEX DOMAINS 



5 



1.1. Light propagation, heuristics and degenerate oscillatory integrals. In [9] the 

second author sketched the main steps of a proof of fll.Sp and gave a full description of 
the geometry of the wave front set. In this work, we provide a complete construction of a 
suitable parametrix for the wave equation, which we then utilize to obtain decay estimates 
by (degenerate) stationary phases. 

Recall that, at time t > 0, one expects the wave propagating from the source of light 
to be highly concentrated around the sphere of radius t. For a variable coefficients metric, 
one can make good of this heuristic as long as two different light rays emanating from the 
source do not cross: in other words, as long as t is smaller than the injectivity radius. One 
may then construct parametrices using oscillatory integrals, where the phase encodes the 
geometry of the wave front. 

In our situation, the geometry of the wave front becomes singular in arbitrarily small 
times, depending on the frequency and the distance to the boundary of the source. In 
fact, a caustic appears right between the first and the second refiexion of the wave front, 
as illustrated on Figure 11.31 and Figure 11.41 (which is a zoomed version at the relevant time 
scale). Therefore, we are to investigate concentration phenomena ("caustics") that may 
occur near the boundary. Geometrically, caustics are defined as envelopes of light rays 
coming from our source of light. Each ray is tangent to the caustic at a given point. If one 
assigns a direction on the caustic, it induces a direction on each ray. Each point outside 
the caustic (and in the sunny side of the caustic) lies on a ray which has left the caustic 
and also lies on a ray approaching the caustic. Each curve of constant phase has a cusp 
where it meets the caustic. 

At the caustic point we expect light to be singularly intense. Analytically, caustics can 
be characterized as points were usual bounds on oscillatory integrals are no longer valid. 
Oscillatory integrals with caustics have enjoyed much attention: their asymptotic behavior 
is known to be driven by the number and the order of those of their critical points which 
are real. Let us consider an oscillatory integral 

(1.8) u,iz) = j^;^ J^ei''^^^^^giz,Ch)dC z eR', C e M, he (0,1]. 

We assume that $ is smooth and that g{-,h) is compactly supported in z and in (. If 
there are no critical points of the map ( — ?■ ^{z,(), so that ^ everywhere in an 
open neighborhood of the support of g{., h), then repeated integration by parts (e.g. non 
stationary phase) yields that |m/i(-2)| = 0{h^), for any > 0. 

If there are non-degenerate critical points, where d(^^ = but det(9^^,^^<l') ^ 0, then 
the method of stationary phase applies and yields ||u/i(z)||loo = 0(1). The corresponding 
canonical form is a Gaussian phase. 

If there are degenerate critical points, we define them to be caustics, as ||u/i(2;) ||loo is no 
longer uniformly bounded. The order of a caustic k is defined as the infimum of k' so that 

\\Uh{z)\\Loo=0{h-^'). 

The most simple degenerate phase beyond the Gaussian is ^f{z, C) = ^ + ^iC + ^25 which 
corresponds to a fold with order k = ^. A typical example is the Airy function. The caustic 



6 OANA IVANOVICI, GILLES LEBEAU, AND FABRICE PLANCHON 



1 

: z3 




Figure 1.2. The caustic for the swallowtail catastrophe. 

is given by 2:1 = and the illuminated side is 2:1 < 0. The next canonical form is given by 
a phase function which is a polynomial of degree 4, namely ^c{z, C) = ^ + + ^2C + ^3 
whose order is k = |; its associated integral is called Pearcey's function and it produces a 
cusp singularity on the caustic which is parametrized by zi = — 3s^, Z2 = 2s^. 

Finally, we conclude this brief overview with the swallowtail integral (which is an os- 
cillatory integral with four coalescing saddle points) whose canonical form is given by a 

polynomial of degree 5, $5(2;, () = y + 2:1 y + 2^2 y + ^sC + ^4- the caustic surface of the 
swallowtail is defined by the condition that two or more real saddle points are equal: it 
is pictured on Figure II. 2[ In the event that two simple saddle points undergo confluence 
when z zq, then the uniform asymptotic behavior of (11.81) contains terms involving the 
Airy function and its derivatives multiplied by powers of /i~2 + 3; the caustic surface is 
smooth [zi < on the figure). If three simple saddles coalesce as z ^ Zq, then the uniform 
asymptotic behavior of (ll.Sp can be described by terms containing the Pearcey function 
and its first-order derivatives, each multiplied by a power of h~^~^i; the caustic surface has 
cusps (two of them in the zi > region on the figure). The swallowtail enters in picture 
when four simple saddle points of (II. 8p undergo confluence as 2; — )■ (which is 2:9 = on 
the figure). 

Such integrals will play a crucial role in the proof of Theorem 11.31 Between two consecu- 
tive reflections of the wave propagating along the boundary, we shall construct a parametrix 
of the form 

uiz, h) = ^ [ et*(^'C'^)^^(z, r/, (, h) dCdrj, 

where the phase is essentially (^{z.C^.rf) ^ ri^ci^X)^ with z = (t,x,y), rj/h is the Fourier 
variable associated to the tangential variable y and C = ^/f] where ^/h is the Fourier 
variable associated to the normal variable x. Note we may restrict to 77 G (1/2,2) which 
corresponds precisely to waves propagating along the boundary and explains the (77, () 
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parametrization for the oscillatory integral. For a particular value zs oi z = {t,x,y), this 
phase will have a saddle point of order 4; it corresponds to 9^$ = = ^ci^sX) ^"^^ 
d(^^c = d'^^c = d^^c = 0: the geometric picture is that of a swallowtail singularity, but 
the decay loss is that of the Pearcey's integral, i.e. h^^'^. For z ^ zs-, our oscillatory integral 
will have only critical points of order at most 3, corresponding e.g. to 9,,$ = = ^c{,zs^ C) 
and di^^c = d'^^c = 0: the picture is, at worst, that of cusps and the loss is that of the 
Airy function, i.e. h^^^. Finally, we notice Figures [L2] and [L^ picture the same singularity 
formation: in 11.21 up to translations, Zi = t, Z2 = —x and z^ = y; Zi < corresponds to 
the (smooth) refocusing wave front in the left part of 1 1.41 while two cusps form on the right 
part after the swallowtail singularity. 

I. 2. An outline of the proof. Let us mention the main ideas of the proof of Theorem 

II. 31 First, we may reduce to the two dimensional the tangential directions will 
produce the usual decay factor when we integrate them out, see Section HI 

Let h G (0,1] be a small parameter {1/h will later be the spectral frequency) and 
< a <^ 1 the distance of the source to the boundary. We assume a to be small as we are 
interested in highly reflected waves, which we do not observe if the waves do not have time 
to reach the boundary. 

From the spectral analysis which will be recalled in Section 13.11 we have an explicit 
representation for the Green function associated to the half-wave initial value problem 
with a Dirac at (a, b) as initial condition at time s: 

(1.9) G{{x, y, t), (a, b, s)) = J2 [ e±*(*-^)v^e^(?^-%efc(a;, r/)efc(a, rj) dr, 

k>i -^^ 

where \k{f]) = ff + V^^^i^k, with —cok a zero of the Airy function and the ek{x,r]) are 
explicit, real- valued functions which are defined in Section 13.11 We now record several 
remarks that will be of help later and relate to various phase space localizations. 

Remark 1.9. We may perform a spectral localization at Xk{v) ~ h~^, which corresponds to 
inserting a smooth, compactly supported away from zero ip2{h\/ Xk{v))'i the flow, this is 
nothing but il)2{hDt) and this smoothes out the Green function. Then we are dealing with 
a semi-classical boundary value problem with small parameter h. With the usual notations 
r = hdt, r] = hdy, ^ = ^d^, the characteristic set of our operator is given by 

= + (1 + x)r/2 

The hyperbolic (resp. elliptic) subset of the cotangent bundle of the boundary x = is 
|r| > |?7|, (resp. |r| < |?7|) and the ghding subset is |r| = |?7|. Fromr^ = {hD^y = h?\k{Dy), 
one gets at the symbolic level on the micro-support of any gallery mode associated to ujk 

(1.10) r/^/3^2/3^^=e' + xr/2. 

Remark 1.10. We may also localize with ipi{hDy), with ipi G C^(]0, oo[), which correspond 
to a Fourier localization along the tangential (e.g. y) direction (notice such a truncation 
is easily seen to commute with the equation, hence the flow). Since we are not interested 
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with waves transverse to the boundary, we may and will assume that on the support of 
4'i{v)4'2ih Xk{ri)) one has k < eh^^ with e small. This is compatible with f ll.lOp since 
cufc ~ /c^/^ and k < eh^^ is equivalent to |^| < e^^^. This fact will later have its importance 
when a < h^f^. 

Remark 1.11. Irrespective of the position of a relative to h, the remaining part of the Green 
function, will be essentially transverse and see at most one reflexion for t G [0, Tq], with Tq 
small (depending on the above choice of e). Hence, it can be dealt with as in [1] to get the 
free space decay and we will ignore it in the upcoming analysis. 

Remark 1.12. Finally, the symmetry of G (or its suitable spectral truncations) with respect 
to X and a will be of great importance: it allows us to restrict the computation of the L°° 
norm to the region < x < a. 

Now, we consider initial data uo(x,y) = ilJ2{h^J — l\g)ilJi{hDy)5x=a,y=o where the ij^j are 
those of remark n.lOi We will use different arguments depending on the respective position 
of a and h. 

The first case is a ^ /i^/'': there, we follow ideas of [6] and write a parametrix for the 
wave equation as a superposition of localized waves for which we can compute the wave 
front set and hence the singularities that appear at different times and locations. The 
construction of [6] has to be significantly altered to allow for the range /i^^^ <C a < /i^^^, 
with a phase which is less explicit but prevents amplifying factors at each reflexion that 
induced the a > h^^"^ restriction in [6]. 

The second case corresponds to data for which the distance a to the boundary is such 
that < a < /i^/^: we write the contribution of our data which is localized in a h^^^ 
cone of tangential directions as the L^(f2) orthogonal sum of whispering gallery modes 
and prove that after a time t the corresponding wave remains frequency localized in the 
same cone of directions of size h^^^, at least up to smooth remainders. While not quite as 
strong as a microlocal propagation of singularities result, this allows for the use of Sobolev 
embedding theorem to recover the "dispersion" by using the size of the Fourier support. 
The contribution of data corresponding to directions with angles with the boundary greater 
than /i^/^ may be dealt with separately, using a crude parametrix construction, as they 
involve only cusp-type singularities. 

Notice that there is an overlap between the two regions: in fact the parametrix con- 
struction obviously provides better bounds in the overlap region, both in size (we gain an 
a^/^ factor in the worst case) and position (the swallowtail occurs exactly once in between 
two consecutive reflexions). Had we reproduced the parametrix construction from [B], we 
would have an epsilon loss in the dispersion estimate because of the a ~ h^^"^ region. We 
thought it was of independent interest to quantify how "far" below h^^"^ the construction 
could be pushed while retaining the most interesting features of p]. 

Remark 1.13. Figure [L3] illustrates the propagation of (part of) the wavefront set of the 
Dirac data; the second picture is a zoomed version of the first one and shows in detail 
the formation of the swallowtail singularity for the part of the wave front moving along 
directions which are initially tangent to the boundary. 
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Figure 1.3. Propagation of the wavefront 



0.08 
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Figure 1.4. The formation of a swallowtail singularity just after the first 
reflection (zoomed image) 

Finally, Theorem II. 31 is obtained for a > in Section [2], Theorem 12. 11 and for a < h^^"^ 

in Section [31 Proposition 13.41 Theorem 1 1.4 1 is obtained in section [21 as a remark at the end 
of the proof of Proposition 12.141 



2. Parametrix for a> hr 

This section is devoted to the construction, modulo 0{h'^), of the Green function in the 
case a > h'^~^. The Green function is represented in proposition 12.41 as a superposition 
of 0(a-i/2) reflected waves. We give a precise analysis of the Lagrangian in the phase 
space associated to each reflected wave. This geometric analysis allows us to track the 
degeneracy of the phases when we apply phase stationary arguments. Our main dispersive 
estimate will be Theorem 12.11 

Let us set h = h/r] and P = {—ihdx)^ + 1 + x — {—ihdt)^. For a > 0, we denote by 
Aq C T*]R the Lagrangian 



(2.1) 



= {(t', r') 36 eR s.t. t' = -2eVl + a + e\ t' = VT 



The set may be parametrized by t'. Let 'ipait') be the unique function such that ipai^) = 

4/9^ (z^ + ^^), from 



and Aa = {{f ji^'ait')}- Let us set p = 1 + a and 6 = then {t 



J\2 
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which we get 



2z' + i = ^i + {t'Y/p' =^ € = p(i +z') = ^{i + vi + {trip') 



and as ip'a = > 0; 

V^: = v^(l+i'V(8p^) + 0(O); 
finally, by integration, as V'a(O) = 0, 

(2.2) Ut') = Vp{t' + + 0{t")) 

2.1. A singular integral representation for the data. We start by a suitable decom- 
position of the smoothed Dirac as an inverse Fourier transform of a superposition of Airy 
functions. 

Lemma 2.1. Let xi £ 6'o, ^o)) with small 6q. There exists a symbol ao{t',h) of 

degree with an asymptotic expansion in h, supported in a neighborhood oft' = with the 
following properties: let 

(2.3) Mt,x,h,h) = ^ f em^-^')+<-+^-e)+^)ei^^^'')a,{t\h)'^^''^''^^ 



2iih J ' ' ' {27rhy 

then Uo is such that, for x > — 1, 

(1) The wave front set of Uq is included in r > 0. In fact, 



WFh{uo) C {r G [Vr+a,ro]}, 
where tq is related only to the size of the support of in t' . Moreover, 

Fuo = 0. 

(2) The initial data uq{0,x, h, K) is a smoothed out Dirac, that is 



(27r/i)2 J 

Proof. Consider the time Fourier transform of Mo, 
Uq{t /h,x,h,h) = I e~^^'^^''uo{t,x,h,h) dt 



[ A2(ri(x + l-r2))et(^»(*')-"*'Vo(t',n)- 
J 



dt' 



Therefore, mq is an average (with compact support in t') of solutions to the equation 

(-9jV + (l + x-r2)/ = 0. 

From dt'lipait') — Tt') = 0, we get r = ipa{t') and therefore there exists 6 such that r 
y/l + a + W, which proves the claim on WF^Iuq). 



DISPERSION FOR THE WAVE EQUATION INSIDE STRICTLY CONVEX DOMAINS 11 

We proceed with the second part of the statement, regarding the initial data, 

[l-KhYlun J 

Let X, s, C) = s{x + 1 — C^) + — t'C, and denote by the set 

The equations defining C,^ read x + 1 + = and 2sC, + 1' = 0. From the first equation, 
we get C 7^ on (recall x > —1). Now, 

■ 2s -2C' 



Hess, 



-2C -2s 



and det(HesScj^^0) 7^ on C^. Therefore is a smooth manifold. 
Denote by vti the projection from to T*M^., that is 

7ri((t',x, s,C) G C^) = = (x, s). 

and by 712 the projection from to T*Rt/, that is 

n2{{t',x,sX) e C^) = {t',-dt4) = (t'X), 

For r' 7^ 0, we have 

n^\t', t') = (f, x = -1 + t"- t'V(4C'), s = -t'/2C, C = t'). 
Therefore induces a canonical transformation from T*Wt' \ {r' = 0} to T*'Rx defined by 
X{t\ r') = {x = -l + - ^'V4r'^ ^ = -t'/2T'). 

Notice that 

x{Aa) = {x = a,^ = e)=T:^,, 

and X is a symplectic isomorphism from a neighborhood of {f, t') = (0, a/1 + a) onto a 
neighborhood of (x,.^) = (a,0). 

The remaining part of the argument is standard: denote by G{t',x) = (f){t',x',Sc,Cc) 
where (sc, Cc > 0) is the unique solution of 1 + x + = Cc t' + 2scCc = 0, then 
G(0, a) = 0, as Sc(0, a) = and Cc(0, a) = a/I + a. By stationary phase in (s, C) we get 

Wo(0,x,/i,n) = ^^-i^ y e*(^(*'^)+'^"(*'))Ao(t',x,;i)ao(t',^)rft' 

where Ao(t', x, /i) is an elliptic symbol of order 0. From dt'G{t', = and ^(0, a) = 

0, we get G(t',a) = —ipait'), and therefore 

x) + ^a{t') = (x - a)Ha{t\ x) with dt'Ha{Q, a) ^ 0, 

and by change of variables = Ha{t' , x) and using that for all F there exists G such that 

we obtain the desired conclusion, since by the above canonical transformation the map 
<7o(t', h) H- G(6, h) is elliptic of degree 0. 
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□ 

Set goit', h) = e:s^"(*')ao(t', h). We proceed with 
Lemma 2.2. Let c > 0, e > 0, then, with p = 1 + a, 
(2.4) sup \go{r/fi,h)\eO{h^). 

Proof. Notice that gQij/K) behaves hke an Airy function from the geometry of so the 
estimate on g^ is really the classical estimate on hi. We provide a direct argument: for 

Ur/m= [ eTM-^'^-'^^ao{t,h)dt, 



and we may integrate by parts using L = (^/'^(t) — r) ^jdt (recall ^a{t) = a/p(1 + t'^/Sp^ 
0{t^))). Notice that 



AT 



«i,Jv(t) 



where aj,jv(t) = t^^-^^i + l3j^pf{t) (by induction, as ^p'^{t) = 0{t) and [N - 2jJ + + 1 > 
[A^ - 2j + 1J+ as well as [[A^ - 2(j - 1)J+ - lj+ > [A^ - 2j + lj+). As such, it remains 
to check that for t G [—1, +1] and a G]0, 1] 

which is trivial if j > N/2 and follows from setting t = ^/as if j < N/2. □ 

Remark 2.3. One may also prove that there exists tq > (related to the support of cq) 
such that 

(2.5) snp\go{T/h,h)\eO{h'^). 

T>TO 

2.2. Digression on Airy functions. We recall a few well-known facts about Airy func- 
tions: let z > 0, 

At{-z) = — [ e'(^'/3— ) ds = A^z) + A_{z), 
27r 

with = A_ and 

-| r+oa 1 3 la 

A_{z) = — / e*(^'/3-s.) _ _L_e-/4e-|i.^ expT{z'/^) = ^e'^'^e--^''^ ^ _{z) 
2vr Jo 2y^2;4 z^ 

with expT(2;^/2) ~ (1 + X]i>iC/^~^) ~ 20r\E'_(z) as z — +oo. We also define '^^{z) = 
^-{z). Therefore, we have 

444 = ^e-t^^'''e^^(^^) with zB = T-t. 

Notice that for u G IR+, B{u) G M and i?(M) ~ J2j>i ^i'""'' +oo. 
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2.3. The parametrix construction. Let -^(C, K) be a function witli compact support in 
C E [1 + c/ii-% Co]- Define 

u{t,x,h) = [ ei(^*+^(^+^-^')+^'/=^)F(C,n)dC- 



One easily cliecks that Fu = and the trace on x = is 

(t,0,n) = ri / ei"^{A+ + A_){h-^/%C^ - l))F{C,h)dC. 



Define / by F{C,h) = J exp{-it'C/h)f{t',h)dt', then 

where J± corresponds to canonical transformations j± on T*]Rf H {r > 1}, 

J±{t\ r ) = it = t'T 2tV^^, t = t'). 
We now set up a few notations: 

• a cut-off function Xoiv) ^ C'(r((l/2, 5/2)) such that Xo = 1 on [1,2]; 

• recaU xi E C^{{-6o,6o)) with smaU 610; 

2 

• let a G [/i3~^,ao], with Oq small; 

• let /5 > be such that a/1 + a — a/1 + af3 > ca, for all a G [0, Oq]; 

• let X2(m) G C°° with X2(w) = for m < /3/2 and X2(m) = 1 for m > /3; 

. let X3(C) e with X3(C) = 1 for 3/4 < C < Co and XsiO = for C > Ci or 
C < 1/2 (with Ci > Co; Co ~ 1 > and small and Ci ~ 1 small). 

Define 

VNit,x,y,h) =-^2^^ J e'T^^UNit,x,h/T])T]Xoiv) dv , 

2'Kh 



UN{t,x,h) =Lj^ f ^jrit<:+six+i-e)+s^/3-m/3ie-i)i+hNB{ie-i)i /h)} 



X X2((C' - l)/a)x3(C)^o(C/^ h) dsdC ., 
v{t,x,y,h)= ^ VNit,x,y,h). 

0<N<Co/^/a 

and let P = - {dl + (1 + x)dl). 

Proposition 2.4. There exists Cq such that the following holds true: 

(1) V is a solution to Pv = for x > —1; 

(2) its trace on the boundary, v(t G [0, l],x = 0) is Oc°°{h°°); 

(3) at time t = 0, we have 

v{0,x,y,h)-{27rh)-' [ e^^^''y+^'-''^^ho{v)Xi{^/v) dvd^ = Oc^{h^). 
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Remark 2.5. Here and thereafter, f{z,h) G Oc°°{h°°) ioi z E T if, uniformly in a G 

Wa,N, 3C„,7V s.t. snp \d^ f {z, h)\ <Ca,Nh^. 

Proof. Obviously v is defined by a finite sum and each is a. solution to Pv^ = 0. We 
postpone the rest of the proof to subsection 12. 5[ □ 

Remark 2.6. We may also define u by a sum from —Cq/ y/a to Co/y/a, and replace t G [0, 1] 
by t G [—1, 1]. The equation enjoys time symmetry and therefore the two points of view 
are equivalent. 

We start by studying u^; from there, we may obtain information on vn by integration 
over rj. This, however, is a non trivial matter, as h = h/t] and integration over rj has an 

effect on exp{iNB{{C - 1)V^))- 
Let 

<j)a,N,n{t,x,t\sX) = {t-t')C+s{x+l-C)+sy3-4N/3{C-lf^+hNB{{C-l)"yh)+Mt'), 
so that 

UN{t,x,h) = \^ j e^^---«X2((C'-l)/a)X3(C)^o(t',^)rft'c?srfC. 

Notice that 

• t' takes values in a compact set close to t' = 0; 

• C takes values in a compact set close to C = 1; 

• the s integral is oscillatory, and as the symbol is independent of s, this yields an Airy 
function (something we will use only to check the trace condition in Proposition 

m- 

Let us set 

Ca,N,h = {{t, X, t', S, C) S.t. dt'(pa,N,h = ds(l)a,Nfi = dc_(t)a,N,h = 0} 

we therefore get a system of three equations defining Ca,N,h, 

C = i^'ait') 

X = - 1 - ^2 

t = t' + 2sc + 4iVC(C' - - Ib'HC - 

Notice that on the support of the symbol in the definition of Mtv, we have C e + (i] 
with Ci ~ L We can thus localize further the symbol with X4:{^) ^ C^, X4 = 1 for 
s G [— Ci) Ci]; as for |s| > \(\, we will have x = — 1 — s"^ < —1 and as such the contribu- 
tion of 1 — X4 will be Oc°°{h°°) (by integration by parts in s) in the x > —1 region. 

Remark 2.7. In fact, one may localize closer to s = 0: if x^i^) = 1 on [—^/Cf — 1, y/(f — 1], 
the same argument provides a remainder term for x > —Eq. Hence, localizing s close to 
implies (i close to 1 and therefore 6q smaller and smaller and the same for oq. 
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We may parametrize Ca^jy^ri by (s, 9) when they are close to the origin: 



X = a + 6*^ - 

t' = -2eVi + a + e^ 

s = s 



C = Vi + a + e^. 



Notice that {s,9) {s,t' = -29VT+a + 9^) is a local diffeomorphism in a neighborhood 
of (0,0), which ensures that Ca^N.h is a smooth 2D manifold. 
Let us denote by Is.a,Nfi the image of Ca^jq^ by the map 

(t, X, t\ S, C) {X, t,^ = d^(j)a,N,h, T = dt(l)a,N,h) 

then Aa,N,h is a Lagrangian submanifold which is parametrized by {s,6): 

t = 2Vl + a + 9^{s -9 + 2NyUT9^{l - ^B'{{a + 6^f^ /K))) 
i = s 

T = Vi + a + e^. 

Lemma 2.8. The Lagrangian submanifold jv.^i is smooth and its parametrization by 
(s, 9) is one to one. 



Proof. One has first to verify that at each point (s, 6), the differential of the map (s, 6) i— t- 
(x, t, ^, r) is injective. But this is obvious since (||, ||) = (1, 0), |^ = and if ^ = 0, 

then ||(s,0) = -2^/TT^. The map is clearly one to one as t{s, 9) ^ t{s, —9) for 9 ^ {]. □ 



We digress for a while and explain how to add the y variable. In the definition of 
f Ar(t, X, I/, /i), we have a phase function 

^a,7V,h(t, X, y, t\ S, C, v)=Vy + V(l>a,N,h/r,it, X, t\ S, C), 

from which we get dy'^a,N,h = ^'^id 



d,^a,NM=y+^a{t')+at-t')+s{x+i-e)+s'/?>+N{e-if\-A/?,+B\{e~i)h/h)) 
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and the full Lagrangian A^^^v./i C T*M^ is the set of points {x, y, t, C,, r], r) such that there 
exist (s, 6, rj) solution to 

x=a + e'^ - 

y = - ^a(-2^(l + a + ^2)^) - 2s(l + a + + 

- N{a + ^2)5(3 + 2a + 2e^){- - B'{{C^ - l)^/h)) 

3 

t =2Vl + a + e^{s -e + 2N{a + 6^)^! - ^B'{{C - l)h/h)) 
^ =r]s 

f] =7] 

T =r]Vl + a + 6^2 . 

Remark 2.9. Notice that for = 0, having t = in Aq^q./i is equivalent to having s = 6. 
This implies x = a and then ?/ = is a consequence of 

(2.6) + a + ^2)5) = -20(1 + a + 9^) + ^6^ , 

and observe that (ES} holds true since ^^(t') = (1 + a + 6'2)i for t' = -26(1 + a + O"^)^ 
and V'a(O) = 0. Therefore we can explicitly compute ipait'), a.s 9 = — + a + ((1 + a)^ + 

2.4. A suitable change of coordinates. We now perform a rescaling of our coordinates 
that provides some useful reductions. 
Set 

t = y/aT, X = aX, y = — t-y/l + a + a^Y, 
and define 'Jn(T, X, Y, h), wn{T, X, h) as 

fAr(t, X, y, h) = 77v(t/A/a, x/a,y + ta~'^l'^\J\ 4- a, h) 

UN{t,x, h) = a~2eh*^^^^WN{t/\/a,x/a, K). 

We define 

Notice that ipa is C°° in (T',a), with support in A/a|T'| < 1 and 

(2.7) ^„(r') = ^^(1 + 0(ar'2)) (recall p = 1 + a) . 

24p2 

Set 

— 1 = az, s = \faa^ t' = y/aT' 

and 

(2.8) C-Vp = a-fa{z) = a 
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SO that 

(t _ t')C + = t^ + ai((T - T')7a(2) + MT')) 

and therefore 



3 

')a,N,h = ty/p + a^(Pa,N,\ 



with 



(2.9) <^a,7v,A(T, X, r, a, z) = ^a{z) (T - TO + i,a{T') + a{X - z) + 

+ iV(-l^l + l5(Azi)), 

3 

where \ = /h will be our large parameter. One may remark that (pa,N,\ is C°° in a and 

(2.10) ^o,^,A = ^(T -T') + ^ + aiX - z) + ^a' + N{-^zl + ^i?(A^i)) 

and we have z > P/2 > on the support of the symbols in our integrals. 
We have now 

(2.11) jN{T,X,Y,h) = j e''^^''wN{T,X,h)r]Xo{v)dv, 
and 

(2.12) WMiT,X,h) = ^^^ [ e'^^--^^x2{z)^^^^XA{V^(T)aoiV^T',h)dT'dadz 

2iT J 2y/l + az 

where we used dt'dsdC, = ^J^^^^ dT'dadz and a?/{2T[h) = y/a\/{2n). Finally, we set 

~ 3 

9 = y/Efx and A = \/r] = a^/h. By our change of variables, the differential operator P 
becomes 

P = a-^Qa with Qa = -d\ - (X - l)dl + 2^1 + adrdy + a9|, 

and 

Q,(e^^^/(T,X)) = e^^^A2Q,/ 

with 

Qa = i^dxY + (X - 1) - 2y^^dT - ai^drf . 
iX iX iX 

Our initial data at T = is now 



j e^4^(''^-*-(^-^)=)xo(r?)xi(v^S/r/) dr/dS , 



(27r/i)2 

and is concentrated at y = 0,X = 1. The new operator Qa has symbol 

a{Qa) = h2 + X - 1 - 2^r - ar^ 
and the positive root in r of cr(Qa) at X = 1 is 

^+ ^Jp + a'E? 



18 OANA IVANOVICI, GILLES LEBEAU, AND FABRICE PLANCHON 

Notice that, as r is preserved by the flow, the bouncing angles at X = are such that 
'^bounce = 1 + 2y^r + ttT^ = 1 + Sq > 1; we are now facing only transverse reflexions, 
however we aim at studying the flow for very large times. 

Remark 2.10. Assuming the worst terms are |^| < y/a, this translates into |E!| < 1 which 
implies r bounded, and for a small Qa degenerates to a Schrodinger operator. 

From Ot' = ^/adf, dg = \fada and 9^ = |(1 + az)"^(9^, we have 

1 - z 



(2.13) aT'<^a,iV,A = ^1(T'^ 



az 



2 



(2.14) d,^a,N^ = ^- z^a 

(2.15) dz^paNX = / {T-T' - 2aVl + az - ANy/^^Jl + az{l - -B'{\z^). 

2\/l + az 4 

The Lagrangian of ifja is parametrized by 



(2.16) T' = -2^^^l + a + a^^\ ^'^{T') = ^ 

^/p+ y 1 + a + a/i^ 

and as 1 + a^; = (^^, and C, = {1 + a + on Ca,N,h, we have 

1 + fi^ = z on Ca,N,h ■ 
In our new set of coordinates, the projection of Aa^N,h onto is, with z = 1 + /i^, 

(2.17) X=l + ^2_^2 

(2.18) Y = 2/i2(/i - a)Hi{a, /i) + '^{a^ - fi^) + 4iV(l - ^S'(A^i))i72(a, /u) 

(2.19) T = 2^p + a/i2(a - /i + 2iVVl + - -B'{Xz'^)) 
where Hi 2 are defined as 



(2.20) Hiia^p)- + ^ 



y/p + ^/p~+af? 



(2.21) iJ2(a,/i) = Vl + /i 



2 3' 9 ' V 3 ' 9^ 



p^y p + ap^ + 1 + |a(l + /i^ 



Remark 2.11. Notice that the parameters are p, a and rj through the A factor in the X, Y, T 
parametrization of Aa^N,h- 

2.5. Proof of Proposition 12.41 We already dealt with the first item. We now address 
the remaining two. 
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2.5.1. The boundary condition. Set 

a 

and recall 

0<7V<Co/v^ 

recall as well that we constructed Fn so that 

which allows to cancel all middle terms in the sum to get 
Let us define 

Ioit,h) = I e^^'-''^U^ih-'/\e-l))e'^^^^'\2X3CTodt'dC 

Hence, it is enough to prove that 

• Iq E Oc°°{h°°) for t>0, uniformly in a; 

• -^A'max ^ Oc^{h°°) for t < 1, uniformly in a, 

Start with Iq. Using our change of scales, /o(t, h) = Jo{a''2t, A) and 



2a/1 + az 

where uiq is a symbol of order 0. As we have dt = a^^/'^dr and X = a^/h > rjh"^^ and 

2 

a > /i3~^, we are left to prove the following: 

Jo(T, A) e Oc7oo(A"°°) for T > , uniformly in a. 

We already computed the derivatives of the phase of Jq. Recall that T' and fj, are related 
by T' = — 2yUA/l + a + a/i^, and a/i^ is bounded. 

(phase of Jq) ^ ~^ 



^/l + az + ^yl + a + aii^ 



„ , , , , , T-T' + 2VivTT^ 

(phase of Jo) = „ ,^ 

2vl + cbz 

where for the first derivative one uses dr' = \fa,dt> and the identity 



7a(^)(T - T') + ^,{T') + -z-2 = a^mt - t'K + 3(C' - 1)^ + Mt') - tVTT^) 
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The first derivative vanishes if z = 1 + fi"^, and the second one vanishes if 

T = -2v/l + a + a/i2(^ + ^1 + ^2) < g. 

As such, the phase has no critical points for T > 0. One has to be careful as the domain of 
integration of the (T', z) variables is very large with small a, as it is like (— c/ ^/a, cj ^Ja) x 
(l/2,ei/a). 

We turn to the details: for z < Zi^, z is bounded. For large |T'|, we get (phase) | ^ 
/i^ T'2, therefore by integration by parts in T' we get decay. If \T'\ is bounded there is 
no critical points in {z,T'). 

We are left with z > zq, where zq is large. We first perform the integration in T'. There 
will be two critical points in T', given by fi± = ±-\/ z — 1. Denote by Ja,±{z) the critical 
values of the phase — 7a(2;)T' + ipaiT'). We have 

Ja,±{^) = (-7a(2)r' + 4(T'))|T'=-2(±)vCTyT+^ 

and 



We are left with the z integral, 

+00 3 

M,)T+iz^+Ja.M^))g^^^^X) dZ, 



where Qa is a symbol of order —1/4, uniformly in a: \dl{z^ga{z, X))\ < Ciz~^ with Ci 
independent of a, 2 > Zq, and Qa is supported in [zq, c/a] with small c (notice we used that 
the mo and Xs terms in Jq are symbols or order 0, uniformly in a). We have 



a,(7a(^)T + -z-2 + J,,±(z)) = ^z± + Ti^{z), 

and for the + case we may integrate by parts in z without difficulties. For the — case, set 

2 T 

J = d,i-fa{z)T + -^t + Ja,-(^)) = - + 



3 V V 2Vl + a2' 

For T > 0, z G C with |Imz| < (5|Rez| and zq <Kez < c/a, we have 

\J\ > Re J > ^ 

'z 



with a constant C which does not depend on a or T > 0. Hence by the Cauchy formula, 
\diJ-^\ < Qz^-^ and J-^ is a symbol or order 1/2, uniformly in a, T > 0. We may then 
conclude by integration by parts in z with the operator X~^dz{J^^^-) (if g is a symbol of 
order m, dz{J~^g) is a symbol of order m — |). 

The remaining integral In^^^^ may be dealt with in a similar way. In fact, the situation 
is easier: on h.a^N,h fl {a; = 0} we have s = ±\/a + 9"^ and therefore 



t = 2Vl + aT¥{±Va + e^ -e + 2A^^axVa + ^(l - -B')) > 3Co. 
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2.5.2. The initial data. Taking into account Lemmata 12. II and 12.21 we are left to prove that 
V]\f{0,x,y,h) G Oc°°{h°°) uniformly in 1 < < Co/y/a for x > 0. Recall x = y/aX and 
from f l2.17p we have 

T = a = /i - 2N^1 + fi^a, with a = 1 - ^B'{\zi) = 1 + C'(^). 

from which we get 

X = 1 + /i^ - = 1 - 4iV(l + fi^)a{N{l + fi^)a - fi) 
and, as iV > 1 and fi^{a^ - 1) = i2^0{{Xzly) = 0((^ - l)/iyz^)) E OiX^^), 

ANy/l + fi^{Ny/l + fi^a - /i) > 2(1 - 0(^)) . 

For A > Ao, Ao large, we get, uniformly in A^, X < —1/2 on the projection of Aa^^^h, which 
is what we need, wn{0, X, h) G Oc°°{h°°) for X > 0, uniformly in A^. We turn to the 
details. As before, we will proceed by integration by parts. 

WNiO,X,h) = — [ e'^^gdT'dadz, 
27r J 

where 5^ is a symbol in a,T', z and 



ij = -ia{z)T' + ^,(T') + a{X -z)+ ay 3 + N{-^z'^ + 



3 \ ' 

For z < Zq, we may localize a to a compact region as d^jip = X — z + a"^ > — Zq and 
large \T'\ will not be a problem. Then for T',a, z in a. compact set we get decay from the 
geometrical observation on the Lagrangian. 

For z > zq, we may again eliminate T' and obtain two contributions. 



iXtp± 



g± dadz, 



where 



2, ,3 , , 4 3 B(Xzl' 



V'i = ±-(^ - 1)2 + o[X - ^) + — + A^( — z^ 



'3^ ' ' '3 '3 A ' 

By integration in a, the case z <X — \ provide decay, while for z > X + 1 we have again 

3 

two contributions ±2/3(z — X)2. The associated phases in z are 

±-{z - 1)2 ± -{z - X)2 + N{--z-^ + ^Y^. 

for which we readily observe that they are non stationary: not only derivatives never 
vanish, but they increase in value with A^ (for A^ = 1 and the two plus signs, we deal with 
large values of z as we have done just above, for the boundary condition). 

We are left with the contributions ofX — 1 < 2;<X + 1, for which again we may reduce 
to compact a as dfji\}± = a"^ + X — z, and we conclude by the geometric observation on the 
Lagrangian. 
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2.6. Decay for the parametrix. This section is devoted to the proof of the foUowing 
result. 

Theorem 2.1. Let a < 4/7. There exists C such that for all h g]0, /iq], all a G [/i",ao], 
all X G [0, 1] , all T g]0, a^^^^] and all Y G M, the following holds true 

(2.22) I lN{T,X,YM<C{2Txh)-\{-^f/' + a'l'h'l') 

0<N<Co/y^ 

Remark 2.12. Note that, in the given range of parameters a,h, the above theorem im- 
mediately implies our main result, Theorem ll.3[ after undoing the rescaling from Section 

El 

We first observe that ■jq(T, X,Y, h) = vo(t,x,y,h) where Vq is a solution of Pvq = in 
X > — 1 with WFhVo C {r > 0}. By Proposition 12.41 (and its proof), the associated data 
at time t = is a smoothed out Dirac at x = a,y = 0. Thus Vq satisfies the classical 
dispersive estimate for the wave equation in two space dimensions and since t = a^^'^T, 
this implies 

|7o(T,X,y,/.)|<C(27r/.)-2(-^)V2. 

Thus we may assume in the proof of (12.221) that the summation is taken over 1 < < 
Recall 

(2.23) MT, X, Y, h) = j e^'^^'^WNiT, X, h)7^xo{v)dv 

where the wjy are defined by 
(2.24) 

27r J 2\/l + az 

We split each in two pieces, Wp^ = Wjy^ + Wn^2', Wn,2 is defined by introducing an extra 
cutoff Xb{.z) G C^(]0, ^o[) in the integral (I2.24p . with zq > 1, close to 1, and Xsl-^) = 1 on 
[P/2, (1 + zq)/2]. wn,i is then defined by introducing the cutoff 1 — x^i^:) in the integral 
fl2.24p . We denote by 'Jn = 1n,i + 1n,2 the corresponding splitting using formula fl2.23p . 
The following propositions obviously imply theorem 12.11 

Proposition 2.13. There exists C such that for all h g]0,/io], all a G [/i",ao]j all X G 
[0, 1], all T g]0, a"-*^/^] and all F G M, the following holds true 

(2.25) I J2 lNA{T,X,Y,h)\<C{27Thy^h^/\ 

2<N<Cq/^ 

Proposition 2.14. There exists C such that for all h g]0,/io], all a G [/i",ao]j all X G 
[0, 1], all T g]0, a'^f"^] and all F G M, the following holds true 

(2.26) I ^ 7^,2(T,X,F,/i)| < C(27r/i)"V/*^/i^/^ 

2<7V<Co/v^ 
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Proposition 2.15. There exists C such that for all h g]0,/io], all ct G [h°',o,o], all X G 
[0, 1], all T g]0, a~^^^] and all Y eM., the following holds true 

(2.27) \^,{T,X,Y,h)\<C{2nhr\{-^y/' + a'/'h'/^). 

The remaining part of this section is devoted to the proof of these 3 propositions. We 
also include a separate section which contains useful geometric estimates, and a section 
where we recall useful known estimates on phase integrals. 

Remark 2.16. The hypothesis a G [h°', ao] with a < 4/7 in theorem 12.11 will be used to see 
that only a few overlap with each others. We will address estimate fl2.22p in the full 
range a G Co] in a forthcoming paper. 

2.6.1. Geometric estimates. We denote in this section by f{a,a^'^) various analytic func- 
tions defined for a and a/i^ small, with /(a, 6) G M for (a, h) G M^. Recall that the projection 
of A.a,N,h onto is given by 

X = 1 + /i^ - 

(2.28) ^ = - <y)Hi + - fi') + 4iV(l - ^B'{Xzl))H2 

T = 2^p + aiJ?{a - /i + 2X^177^(1 - ^^'(A^i)) 

with 2 = 1 + /i^ and Hi^2 of the form (see fl2.20p ) 

ffi = /o(a,a/i2), /o(0,0) = l/2 
^^■^^^ i72(l + /i')^'/' = /i(a,a/i') + /iV2(a,a/i2), /i(0,0) = l/3, /2(0, 0) = -1/6 . 

Let us rewrite the system of equation fl2.28p in the following form 
X = 1 + /i^ - ^2 

(2.30) ^x_o,.2^,. , 2/ 3 _ „3n , 9 7T /I , „2n-1/2 ^ 



Y = 2^2(^ - a)H, + -{a'- ^^) + 2H,{1 + fiT'^'i = - a + /i) 

3 2y/p + ai2^ 



and 



(2.31) 2iV(l - Is'iXz'i)) = (1 + p'r'/'i ^ -a + p) 

Then (12.301) and (I2.3ip is obviously equivalent to (12.280 . For a given a and a given point 
{X,Y,T) G M"^, (I2.30p is a system of two equations for the unknown {p,cr), and we will 
use the fact that (I2.3ip gives an equation for N. Recall that we are looking at solutions of 
(I2.30p in the range 

a G [/i", ao], a < 4/7, a\p\^ < eo, < T < a'^^^, X G [0, 1] 
with oo, eo small. Let us denote i? = 2(1 — 3Y/T). 
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Lemma 2.17. Let T > Tq > 0, X e [-2,2], F G M . There exists i^j{X,Y,T,a) G C, 
j = 1, 2, 3, 4 such that 

{fi G C, a|/ip < eo, 3cr G C, {fi, cr) is a solution of l{2.3(]\) } C {fii, /i2, /is, /i4} 

Moreover, there exists a function f*{a, afi^) with /*(0, 0) = 1, and constants Co, Ci, C2 > 0, 
Rq, Mo > such that the following holds true: 

(a) If \R\ > Rq, two of the fijf^:{a, a/i^) are in the complex disk D{\/R, A), the two others in 
the complex disk D{—\fR, A) with A = Co(l/T+ "^^1-^^-* ). Moreover, one has \/\R\ > 2A. 



(h) If\R\ < -Ro (^''^d \R\T > Mq, two of the fijf^,{a,afij) are in the complex disk D(\/R, A), 
the two others in the complex disk D{—^/R, A) with A = . Moreover, one has \/\R\ > 

2 A. 

(c) If \R\ < Ro and \R\T < Mq, one has < C^T'^I'^ for all j . 

Proof. We first get rid of a. The second equation in (12.301) is of the form Y = Bq + aBi + 
cr^B2, thus by the first equation, we get 

(2.32) Y-Bo = (r{Bi + B2il+fi^-X)) 
and then the first line of ( ]2.30p gives an equation for fi 

(2.33) {Y - B^f = (1 + ^2 _ ^^(^^ ^ ^^(^ ^ ^2 _ ^))2 ^ 

where 

B, = 2fi^H^ - 2/^73 + 2H2{1 + fiY'^\^L== + fi) 

2y'p + an'' 

^^■^^^ B^ = -2^,'H^-2H2{l + ^^'r'/' 

B2 = 2/3 . 

By (12.201) and (12.291) one gets through explicit computation the identity 

(2.35) yw, fo{0,w) + f2{0,w) = l/3 

and this implies, (we use aiJ^^^ = {afi'^)fi'') 

Bo = -/i'r/6/3 + 2/3/X/4 + T/3/5 

5i + ^2(1 + /i^ - X) = -2X/3 + /e 

with /z(0,0) = 1 for / = 3,4,5 and /6(0,0) = 0. Let D = -2X/3 + /e. We may then 
rewrite (12.331) as 

(2.37) P+P_ = 36(1 -X)^ 

with 



(2.38) P± = {f.fxf - ^t^(2/4 T m/f* -R', R' = R+ /r, AlO, 0) = 

and/, = y73,/,(0,0) = l. 
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By classical arguments on perturbations of polynomial equations, fl2.37p implies that for 
a and a|/ip small, the /i equation fl2.33p admits at most 4 complex solutions (at most since 
we have the constraint a|/ip small). Set /i* = yu/*(a, a/i^). Then yU, H- /i* is a holomorphic 
change of coordinates and \^ — 1| < Cte(ao + eo). With the notation g± = (2/4=f3Z})//^,, 
the two roots //^ ^ of satisfy the equation (recall that and R' are functions of (a, a/x^) 



(2.39) fxl^, = gjT±y/R' + igjTy, e = ± 

Assume first \R\ > Rq with _Ro large. Then by fl2.39p and fj^O, 0) = 0, there exists Co such 
that 

(2.40) /i;.,GD(v^,Ao/2), /il^,GD(-v^,Ao/2), = Co(l/T + ^^i±M) . 

For Y^li?! > 2Aq, the two disks D± = D{±\/R, A^) do not overlap and dist{D^, D_^) > 
\/\R\- Thus there exists C3 such that 

\P+P^{fi)\>C3\R\Al>C3C^,\R\/T^, V/i* gC\(D+UZ}_), 

and this contradicts f l2.37p for > i?o large enough, and proves (a). 

For \R\ < Rq, and |-R|r > Mq, with Mq large, f l2.40p remains true and thus we get 

(2.41) D(v^,Ai/2), ^,*_^^eD{-^/R,A,/2), A, = Co{l/T + ""^^ ) . 



Set A2 = — H=. Since \R\ < Rn and T < a~^/^, for C large enough one has Ao > 2A\. 

tJ\R\ I I _ u _ , to to / _ i 



The two disks D± = D{±Vr, A2) do not overlap and dist{D+, D_) > ^/\R\ for Mq > 2C. 
Thus one has 

\P+P-if^)\ > C^AllRl = CiC^T^, V/i* G C \ {D+ U D^) , 

and this contradicts fl2.37p for C large enough, and proves (b) for Mq large enough. Finally, 
for \R\ < Ro, and \R\T < Mq, one has clearly by (ESU) and T < a'^/^ < CgT-^/^ 

and thus for c > 2C5 

V/i G C such that > cT~^^\ |^+^-(/i)| > Cec^T^ 

This contradicts 02.370 for c large enough. The proof of Lemma 12.171 is complete. □ 

Let us now study the equation fl2.3ip . which provides . Since B'{u) G 0{u~^), 
z > (3/2> 0, N < Cofl-i/^ and a > h'^'\ one has 

(2.42) \NB\\z^\ G 0(iVA~2) = O^Nh^/a"^) G 0{a~'"^h'') G 0(1) . 

Let < /i >= (l + /i2)V2. From(T2_^2 ^ l-X, we get for X G [-2,2] |cr-/i|/< /i > G 0(1). 
Therefore ( I2.3ip implies 

(2.43) 2A^ = r<l>,(/i) + 0(l), $,(/i) ^ 



2 < fi > a/p + a/i^ 
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Let [/ = {/i G C, < 0.5 or |/m(yu)| < |i?e(/i)|/v^}. Then $a(Ai) is bounded on U and 
(2.44) |<|.„(/i) - <|.„(/i')| < ,^A a + „ J , V/i,/i' G ^. 



< >< l/i'l >■ 

Observe that for 6 G M, and |6| > 2r, the complex disk -D(6, r) is contain in [/. 

For a given point (X,F,r) G [-2,2] x M x [CCoo-^/^]^ let us denote by ^^{X,Y,T) 
the set of integers > 1 such that (12.281) admits at least one real solution (/x, a, A) with 
a|/ip < eo and A > Aq . We denote Af^{X,Y,T) the set of complex such that fl2:28ll 
admits at least one complex solution {fi, a) with fi E U and a|/ip < eo and A > Aq. Observe 
that A/'(X, y, T) depends on a . For E G N, \E\ will denote the cardinal of E. Observe 
that (I2.43P implies for an absolute constant Nq 

(2.45) Af{X, Y, T) C [1, T/2 + A^o] ■ 

Lemma 2.18. There exists a constant Cq such that the following holds true. 

(a) Forall{X,Y,T) G [0, 1] xMx [0, fl-^s]^ one has \M{X,Y,T)\ < Co, andM^{X,Y,T) 
is contain in the union of 4 disks of radius Cq. 

(b) For all (X,F,T) G [0, 1] x M x [0,0-^/^]^ the subset ofN, 

J\fi{X,Y,T)= \J Af{X'X,T') 

\Y'~Y\ + \T'-T\<1,\X'-X\<1 

satisfies 

\Af,{X,Y,T)\<Co. 

Proof. We start with (a), which is a consequence of (I2.43p . since by lemma [2. 171 for a given 
(X, Y,T > To), there is at most 4 possible values of /i (for T < To we use (12.451) ). 

We proceed with (b). By (I2.45p . we may assume T > Ti with Ti large. Recall i? = 2(1 — 
3Y/T). Let {X',Y',T') such that \Y' -Y\ + \T' -T\ < 1,X' G [0,1]. Set i?' = 2(l-3r7T'). 
One has |-R — R'\ < C{1 + \R\)/T. Let us first assume > 2i?o, with i?o as in lemma 
[2Tn Since T is large, one has \R'\ > Rq, and \R'\ ~ \R\. Let N' G N{X',Y',T') and /i' 
such that (12.281) holds true. By lemma [27171 (a), one may assume /i'* G D{\fR , A). Take 
/i* G D{\/R,A) associated to (X, F, T). Since is real, one has R' > Rq, hence R > 2Rq, 
and therefore /i G ?7. Let N G A/''^(X, F, T) associated to /i. From a^^"^ < 1/T one gets 

(2.46) l/i - /i'l < - < C(A + A'+\VR- VR\) < '^^^^^^ . 

By (12.431) and (I2.44p . this implies since a\R\ ~ a|/i'^| < eo, 

(2.47) 2\N-N'\ < |T'-T|$,(/i')+T|$,(/i')-$a(/i)|+0(l) 

<C(«+1/I^l)^^ffl + 0(1)G0(1). 

Let us now assume \R\ < 2Rq and T|i?| > Mo + 8. From \RT-R'T'\ = \2{T -T') -6{Y - 
Y')\ < 8, we get \R'\T' > Mq. We may thus apply lemma EH] (b). Let N' G N{X', Y', V) 
and /i' G M such that (I2.28P holds true. Since is real one has R! > 0, thus R'T' > Mq, and 
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this implies R>0 (take Mq large). Moreover one has | | < C(l + |i?|)/T, < 3/2o, 
and also \R'\ ~ \R\. By the same argument as above, we get now |/i — < and 

since + < C^/\R\, one gets l/i^ - < C(l + Ro)/T and thus by (ICTj) and 

(2.48) 2|iV - N'\ < CT{a + 0(1))(1 + Ro)/T G 0(1) . 

Finally, for \R\ < 2Ro and T\R\ < Mq + 8, one has T'\R'\ < Mq + 16, thus by part 
(c) of lemma [2.171 one has < CT^^'^, < CT~^^'^, and thus we get in that case 
l/i^ - < C/T thus (12:481) holds also true in that case. Since G 7V^(X,F,T), ^M), 
fl2.48p . and part (a) of our lemma imply (b). The proof of our lemma is complete. □ 

2.6.2. Phase integrals. We first recall the following lemma, for which we refer to [TB] . 

Lemma 2.19. Let K G M. a compact set, and a{C,,X) a classical symbol of degree in 
A > 1 with a{C,, A) = for ^ ^ K . Let k > 2, cq > and $(0 o, phase function such that 

2<i<fc 

Then, there exists C such that 

I J e^^*(«)a(e,A)| < CA-l/^ VA > 1 

Moreover, the constant C depends only on cq and on a upper hound of a finite number of 
derivatives of^^'^\a in a neighborhood of K. 

The next lemma will be of importance to us. As such, we have included its proof for the 
sake of the reader while not claiming any novelty. 

Let -ff(0 be a smooth function defined in a neighborhood of (0,0) in M^, such that 
H{0) = and Vif(O) = 0. We assume that the Hessian H" satisfies rank(if"(0)) = 1 
and Vdet(iJ")(0) 7^ 0. Then the equation det{H"){C,) = defines a smooth curve C near 
e with G C. Let s — )■ ^(s) be a smooth parametrization of C, with ,^(0) = 0, and 
define the curve X{s) in by 

Lemma 2.20. Let K = G M^, |^| < r} , and a{^,X) a classical symbol of degree in 
A > 1 wtth a{^, X) =0 for l^ K. Set for x G close to 

(2.49) J(x, X) = j e'^("-«-^(«»a(e, A)^^ • 

Then for r > small enough, the following holds true: 

(a) If X'{0) 7^ 0, there exists C such that for all x close to 

|/(x,A)| < CA-^/^ 

(b) If X'{0) = and X"{0) 7^ there exists C such that for all x close to 

|J(x,A)| < CX-^'\ 
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Moreover, if a is elliptic at ^ = , there exists C such that 

|/(0,A)| > C'\~^'\ 

Proof. By a linear change of coordinates in ^, we may assume H{^) = ^i/2 + 0{C,^). Set 
^) = x.^ — H{^). Then <l>^^(x,^) = x — if^^(^). Therefore, there exists a unique non 
degenerate critical point ^i{x,^2) in the variable ^i, and the critical value \l/(x,^2) satisfies 

and by stationary phase in ,^1, one has 

By lemma [2.191 it remains to prove: 

(a') If X'(0) ^ 0, there exists cq > such that for all (x,^) close to (0,0) > cq. 

(b') If X'(0) = and X"{0) ^ 0, there exists Cq > such that for all (x,^2) close to 
(0, 0) I^I^GI > Co, and in the case a elliptic the lower bound at x = 0. 

Let us prove (a'). Since G{x, ^2) is smooth and r small , it is sufficient to prove d'^^G{0, 0) 7^ 
0. The Taylor expansion of H at order 3 reads as follows 

H{0 = ei/2 + + + c^ie2 + ^^2' + 0{^') . 

Thus one has ^i^(0,6) = -c^| + 0(^f) and we get -G(0,6) = 3rf^| + 0(^|). Thus 
(9|2G(0, 0) 7^ is equivalent to d 7^ 0. 

On the other hand, one has detif"(^) = 2cC,i + QdC,2 + 0{^'^), and since by hypothesis 
Vdet{H"){0) ^ 0, one has {c,d) ^ (0,0). Moreover, one has 

X{s) = H\as)) = {Us) + 0{s'),0{s')) 

and therefore ^'(0) 7^ is equivalent to ^i(O) 7^ 0. This in turn is equivalent to the fact 
that ^1 is a parameter on C, which is equivalent to d ^ 0. 

Let us now prove (b'). Since ^'(0) = 0, we get d = and therefore c 7^ 0. Now, ^2 is 
a parameter on C, and we have ^1 G 0(^|) on C. We will use a Taylor expansion of H at 
order 4, but since ^1(0,^2) is quadratic in ^2, and ^i(s) quadratic in s, we will just need 
the ^2 term, i.e 

H{0 = ei/2 + + be^^2 + c^iil + 6^2 + - + o{e) ■ 

Then we get deti/"(0 = 2c^i+4(3e-c2)e| + 0(e^). Therefore ^1 = 2{c-^e/c)il + 0{Q) 
is an equation for C, and we get that X"{Q) 7^ is equivalent to X('(0) 7^ 0. this in turn is 
is equivalent to 7^ 2e. On the other hand, we easily get — G(0,^2) = (4e — 2c^)^| + 0(^|). 
Finally, for a 7^ and 6(^2, A) a symbol of degree elliptic at ^2 = 0, and supported in 
1^2! with r small enough, one has clearly 

I j e'K<t+oimi,{^p^^x)di2\ > C'X-'/' 
which completes the proof. □ 
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2.6.3. Proof of proposition \2.1^ Recall 
(2.50) 

^;v,i(T,X,n) = / e^^^^'-'-X2{z f'y,^^\ ,{V^cT)ao{^T\ h){l-x,){z) dT'dadz 

2n J 2vl + az 

where the phase fa,N,x is defined by (see (12. 9p ) 

(2.51) 'fa,NAT, X, T', a, z) = ^a{z){T - T') + ^,{T') + a{X - z) + a^S 



O A 



For ej = ±, define 



3/2 



(2.52) ^N,e„e,{T, X, z; a, A) = 7a(^)T + 2ei/3(z - 1)=^/^ + 2e2/3{z - X) 

-Ni^l-^BiXzl)). 

Lemma 2.21. The following identity holds true 

(2.53) WN^,{T,X,h) = Y,j e*'*"-i-2e,,,,,(^;a,A)rfz + i?^,„(T,X,n) 

where Be .(z;a, A) are smooth functions of z with support in [(l + 2;o)/2, (C^ — l)/a] (remark 
that C,i> 1 is a upper hound for the support of Xs)- Moreover, 

\z''d[Qepsiionj \ < Ciz"^^"^ with Ci independent of a,X. 

The remainder Rj^^a{T, X, h) is Oc°°{K^) for X E [0, 1], T E [0, a^^^^], uniformly in a, N. 

Proof. The proof is a simple application of stationary phase in (T', a) in the integral (12.501) . 
Recall z > {1 + Zo)/2 > 1 > X on the support of (1 — X5){z), and z < {(f — l)/a on the 
support of X3(a/1 + az). The a integral is equal to 

(2.54) Ji = j e'^("'/3-<x{.-x))^^^^^)^^ 

One has y/a{z — X)^/^ < — 1. Thus, by stationary phase near the two critical points 
s = ±1 and integration by part in s elsewhere, we get 

(2.55) Ji = A-i/2(^ - X)-i/^(e2/3*"(^-^)'''6+ + e-^/^^^^^-^^'^'fe.) + 0(A-°^(^ - X)-°°) 

where b±{y/a{z — X)^/^,A(z — X)'^^'^) are symbols of degree in the (large) parameter 
A(z — X)^/^. Next, the T' integral is equal to 

(2.56) J2 = j e^^(-^"(^)^'+'^''(^'))ao(xAT', h)dT' . 
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Recall that 

H'^ + 1-Z 



a/1 + az + ^1 + a + a/i^ 



Thus we get two distinct critical points = =F2-\/^^^vT+~a^. The associated critical 
values are —'ya{z)T'^+il)a{T'^) = ±2/3(^ — 1)^/^. As before for the a integral, we perform the 
change of variable T' = s\/z — 1, in order to have the two critical points s± = =fV1 + az 
uniformly at finite distance in z. One has y/a{z - 1)^/^ < ^(f - 1, and using ( 12 .yp and 
(12. 8 p we get, again by stationary phase near the two critical points s± and integration by 
part in s elsewhere, 

(2.57) J2 = X-'/\z - l)-i/4(e2/3.A(.-i)3/2^^ ^ g-2/3.A(.-i)3/2^_) _^ 0(A-°^(z - 1)"°°) 

where c±{y/a(z—lY^'^, \(z—lY^'^) are symbols of degree in the large parameter \(z—l)^^'^. 
By (12.550 and (12.570 . one gets that formula (12.531) holds true with symbols 

which completes the proof of Lemma 12.211 □ 



Let 



(2.59) 



In order to prove proposition 12.131 we are reduced to proving the following inequality: 

(2.60) I lN,i,e,,eAT,X,Y,h)\<C{2nh)-^h'/' 

2<N<Co/y/a 

with a constant C independent of h g]0, Hq], a E [/i", ao], X E [0, 1], T E [0, a^^/^]. 
For convenience, we take Z = z^^"^ as a new variable of integration so that 

(2.61) WN,i,er,e2{T,X,h) = j 6^^''^'^^'^^ e,,,,,{Z ; a, \)dZ ; 

Qei,e2{Z; a, A) are now smooth functions of Z with support in [((H-zo)/2)^/^, ((Ci — l)/a)^''^]- 
Since dz = 2Z-^/^dZ/3, we get |Z'9^e| < QZ^^/s ^j^j^ independent of a, A. One has 
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(2.62) 



H,^^^^,^ =Z-^'^{^{1 + aZ2/3)-V2 + ,^(^2/3 _ 1)1/2 ^ ^^^^2/3 _ ^)l/2 

+ ei(Z2/3 - l)-i/2 + e2X(Z2/3 - . 

We will first prove that fl2.60p holds true in the case (61,62) = (+,+)• From fl2.62p . we 
get that the equation dzHa,+,+ {Z) = admits an unique solution Zg = Z^(T,X,a) > 1, 
such that 

lim Z^{T, X,a) =1 uniformly in X, a 

(2.63) > lz'J'dlH^,^AZ,) = - f (1 + aZf r/\l/2 - aZf) 

- 1/2(^2/3 - l)-3/2 _ X/2(Z2/3 - X)-'/^ . 

Therefore, the function Ha^+^+{Z) is strictly increasing on and strictly decreasing 

on oo[. Observe that 

(2.64) H,,+,+a) = |(1 + ay'/^ + (1 - Xy/^ lim H,,+AZ) = 2 • 
For all k one has 

(2.65) sup |9|(iV5'(AZ)| < CkN\-^Z-^''+^^ < C'^h^a-'"'^ < C'^h" , v = 2- la/2 > . 

Z>1 

Let To > 1. We first prove that (^M) holds true for T e [0,To]. Since Ha,+,+ {Z) < 
C{1 + T), for N > N{To) = C{1 + Tq), one gets |(9z$jv,+,+(^)| > cqN with cq > 0, and 
\d^dz^N,+,+ {Z)\ < CkNZ~^ for k > 1. Therefore, by integration by parts in Z in fl2.6ip 
with the operator L{Q) = X~^dz{{dz^N,+,+)~^Q), one gets an extra factor {XNZ)~^ at 
each iteration. Thus, we get W]y^i^+^+ G 0{N~°°X~°°), and this implies 

(2.66) sup I V 77v,i,.„..(T,X,y,/i)| gO(/i~). 
T<n,xmilYm ^^T,)<N<c,/v-a 

Next, for T G [0,To] and 2 < < N{Tq), one may estimate the sum in fl2.60p by the 
sup of each term. But in that case, we know by (12.621) . f l2.63p . and f l2.65p that there exists 
at most a critical point of order 2 near Z = Zq for $Ar^+^+, and 

\dz^N,+,+ \ + |9|$iv,+,+ | + \dl^N,+,+ \ > . 

Moreover, by the second item of (12.640 . and N > 2, one has a positive lower bound for 
\dz^N,+,+ {Z)\ for large values of Z; thus, large values of Z yield 0{X~°°) contributions 
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to m;7v,i,+,+5 and eventually the worst contribution to wn^i^+^+ will be the critical point of 
order 2 near Z = Zq. This provides 

|wjv,i,+,+ (r,X, h)\ < CA"^/^ with C independent of T G [0,To],X e [0, 1] . 

Since ai/2^"i/3 = h^/^ = {h/rifl^, we get 

(2.67) sup I V iNXeueAT.X,Y,h)\<C{To){2nhy''h^'\ 
T<To,xe[o,i],reK 2<iv<^(To) 

Next we prove that that f l2.60p holds true for T G [Tq, a^-^/^]. As before, we may assume 
N < CiT, with Ci large, the contribution of the sum CiT < N < C^j ^Ja being negligible. 
Recall that we have Z > Zq = ({1 + 2;o)/2)^/^ > 1 on the support of €>+,+ in formula 
f l2.6ip . By the first item of fl2.63p . one may choose Tq large enough so that Z^{T, X, a) < 
(1 + Zq)/2 < Zq for all T > Tq. By the last item of fl2.62p . increasing Tq if necessary, and 
using (12. 65 p . we may assume with a constant c > 

(2.68) \dl^N,+,+ {Z)\ > cTZ-^^\ \/Z > Zo, VT > Tq, ViV < Coa-^^^ 

Therefore, on the support of 6+,+, the phase $Af,+,+ admits at most one critical point 
Zc = Zc(T, X, N, A, a) and this critical point is non degenerate. Since N > 2, from the first 
two items of (I2.62p we get zl^'^ < T, and this implies zl^^ ~ T/N. If T/N is bounded, Zc 
is bounded, and since dz^ < — c < for large Z, we get by stationary phase 

\wN,i,+AT,X,h)\ < CX~^^^T~^^^ withC independent of A^. 

li T/N is large, then we perform the change of variable Z = s{T/N)^ in (I2.6ip : the unique 
critical point Sc remains in a fixed compact interval of ]0, oo[, one has < —c(T/N)'^ < 
for s large and also 

a,^e+,+ (s(T/A^)^a,A) < Cfe(iV/T)2s-2/3-^ 
Thus, by stationary phase 

(2.69) VT G [To, , sup sup |m;^,i,+,+(T, X, ^)| < CA'^/^T-^/^ 

2<Af<CiTXe[0,l] 

By lemma 12. 18^ we know that for any given M = {X, Y, T) there is at most Cq values of N 
such that the projection of Aa,N,h intersect the ball of radius 1 centered at M; therefore, 
we will prove that the previous arguments imply 

(2.70) sup I Yl 77v,i,+,+ (r,X,F,/i)| <C(ro)(27r/i)-2(a-i/4/,i/2) 
re[ro,a-i/2l,xe[o,i],yeR 2<n<Co/V^ 

and since a~^^^h^^'^ < h^^^, we will get that (I2.60p holds true. Let us now explain more 
precisely how one can estimate the sum in fl2.70p by the supremum over A^. 
Let Gn{T,X, X,a) = ^n,+,+ (T, X, Zc(T, X, N, X,a)] a, X). The stationary phase at the 
critical point Zc = Zc(T, X,N, X, a) in (I2.6ip gives 

(2.71) WN,i,+AT. X, h) = A-^/2r-^/'e*^^^(^'^'^''^V7v(T, X; a. A) . 
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By fl2.55p . (12.571) . (12.581) . we know that the iPn{T, X; a, A) are symbols of degree in A, and 
d^i^N < Ck\^^ with Ck independent of 2 < iV < CiT. This gives with A = a^''^/h = X/rj 

(2.72) 7JV,i,+,+(T, X, y, h) = 



This is an integral with large parameter A and phase 

L^(T, X, Y, 7], A) = v{Y + Gm{T, X, At^)). 

By construction, the equation 

dr,LN = Y + Gn{T, X, A, a) + A9aGjv(T, X, A, a) = 

implies that {X,Y,T) belongs to the projection of Aq at/i on M.^. Let T G [To,a~^/^],X G 
[0, 1] and F G M be given. For N ^ Afi{X, Y, T), one has therefore d^LN{T', X', Y', r], A) ^ 
for all A and all X' G [0, 1], \Y' — Y\ + \T' — T| < 1. This implies, since dr/L^ is linear in 
Y, |5^Ljv(T, X, y, ?7, A)| > 1. Moreover, one has, with independent of N,T,X,ri, X,a 

(2.73) |aJ(a,L^)l <Cfc. 

To prove (12.731) . we just use that dxZ^ satisfies 

dxZ,dl<l>M,+AZc) = ~dxdz<l>N,+AZc) = NZ,B"{\Z,) 

and thus from fl232D and ([MS]), we get for all A; > 1, {rid,^fZ^ = {XdxfZ^ G Oih"). Then 
(I2.73P follows from 



N 

XdxGNiT,X, A) = A(9A$^,+,+)(T,X,Z,;a, A) = -(-fi(AZ,) + XZ,B'{XZ,)) . 
Therefore, by integration by parts in t] in (12. 72 p . we get 
(2.74) sup I lNx+AT,X,Y,h)\eO{h^). 

TG[To,a-i/2],Xe[0,l],yGR N<^AfjiX,Y,T) 

Finally, by lemma 12.1 8[ one has \N'i{X,Y,T) \ < Cq, and therefore, we get from (I2.74p and 
f l239|) that fl2T0|) holds true. 

Next, we show that (I2.60p holds true for (61,62) = (—,+). In that case, from the last 
item of (I2.62P and X G [0, 1], one gets dzHa-,+ < for T > 0. Therefore the function 

Ha,-,+{Z) decreases on [1, oo[, from i/a,-,+ (l) = + (1 - X) Ha,-,+{^) = 0- 

The equation dz^N-,+ = admits an unique solution Z^, and this critical point is non 
degenerate. We can thus argue as we have done before for the (+, +) case. 

Finally, one sees that the case (61,62) = (+, — ) is similar to the (+, +) case, and the 
case (61, 62) = (— , — ) is similar to the (— , +) case. We leave the details to the reader. 

The proof of Proposition 12. 131 is complete. 
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2.6.4. Proof of proposition \2. 14\ Recall that 



(2,75) u,.,(r. X, n)^tfA[ ^3((1 H- "f) 

27r J 2(l + az)2 



as well as dT"Pa,N,x = , ^ ^r^= ^^ and d„Lpa,N,\ = X - z + a^. 

Let K = {T' = 0, cr G [—1, l],z = 1}, let w be a small neighborhood of K and Xei^', a, z) G 
C^{u) equal to 1 near K. Since for z in the support of the integral (12. 75 p . one has z G 
[(3/2, zq] with zq > 1 close to 1 , decreasing 2;o > 1 if necessary, we get by integration by 
parts in T',a 

(2.76) ' JvTTTT.) 

2v 1 + az 



with 0(A~°°) uniform in T, X, A^, a. Moreover, xiT' , a, a, h) is a classical symbol of degree 
in with support (T', a, 2) G w and a is just a harmless parameter in x- 
We first perform the integration with respect to z in (12.761) . Recall 



iPa,N,x{T, X, r, a, z) =ia{z){T - T) + ^,(T') + a{X - z) + a^S 

1 



(2.77) 



+ N{--z-^ 



X 



B{Xz-2)) 



dz(Pa,N,x{T, X, T',a, z) 



^ ^' --a- 2Nz^'\l - -5'(A^3/2)) 



2(1 + 0^)1/2 



Thus, (pa,N,x admits a unique critical point zdT, T', a, a, A), and we are just interested with 
values of the parameters such that Zc is close to 1. With u = {^^-y- — o-)/2N, this means 
u close to 1. Since cr is close to [—1,1] and > 1, we may thus assume T = T/AN 
close to [1/2,3/2], say T G [1/4,2]. We denote by g(T,T' ,a; a, N, X) various functions 
which are classical symbols of degree in A and with parameters a, A^; in particular, with 
w = (T,T',a), for all a, there exists Ca independent of a,N,X such that \d!(^g\ < Ca for 
all f G [1/4, 2] and all (T', a) close to (0, 0). 

We denote by fk{a,T,T'/N,a/N) functions which are homogeneous of degree k in 
(T'/N,a/N). The notation Oj means any function of the form / = 'Yl,k>j fk- We will 
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use the following functions 

-TO 



1 + + 4af 2 



(2.78) 



1 _ P„ 



VI + a + VI + aFo 



i^i = -^(T72 + ^(l + aFo)i/2). 

Lemma 2.22. (a) One has 

Zc = Fo + Fi + 02 + go/X^ . 
(b) The critical value '^a,N,x{T,X,T',a) = ipa,N,x{T, X,T' , a, Zc) is equal to 
^a,iv,A ={X - Fo)a + + HoT' + tP,{r) 

(2.79) 1 

-Y^(3^72 + ^)' + aiV03 + WA' 

+ iV(4T7a(i^o) - ^i^o^' + ^72(T, a, iV, X)/X^) . 
Proof, (a) The equation for Zc when A = oo is 

^1/2(1 + a^)i/2 = ^ _ ^(T72 + a(l + a2)i/2) ^ 

The solution of this equation is clearly of the form z = fk{o,, T, T'/N, cr/N) with /o = Fq 
solution of Fo(l + aFo) = T^, and we get Fi by a Taylor expansion at order 1. Then (a) is 
a consequence of the implicit function theorem applied to 

.V2(i + _ h'iXz'/^)) =f- ^{T'/2 + .(1 + a.)V2) . 

To prove part (b), one may of course insert the formula for Zc into the definition of ipa,N,x- 
Another way is to use 

dai-^a,N,X -dX - -i'a)=- 

(2.80) ; . ' 

5r'(*a,7V,A -aX - ay 3 - ^a) = - la{Zc) ■ 



Using part (a), this system is easily seen to be integrable and yields formula f l2.79p up to 
an integration constant which is easy to compute when T' = a = 0. Moreover, when a = 
and A = oo, one has zl^'^ = T — i^{T' /2 + a); therefore, one can easily compute the value 
of the critical value when a = 0, A = oo, and this provides the first two terms on the second 
line of (12. 79 p . The proof of our lemma is complete. □ 
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From aT G 0{a^/^) and iVA"^ e 0{h^), we get dl^a,N,\ = -Nz'^''^ + 0{a^''^ + /i^) for 
any z close to 1. Decreasing qq and if necessary, we get by stationary phase, 
(2.81) 

I e'^'^^'^'^xiT' , a, z; a, h) dT'dadz = -^L= [ e*^*'"-^-^ .x{f, T\ a; 1/iV, a, K) dT'da . 
J V AA^ J 

Here, x is a classical symbol of degree in h, with harmless parameters a,l/N. Let us 
define Jn,2(T, X, Y, h) by the following formula, where A = cr'l'^ jh = X/t]: 



(2.82) 7;v,2(T, X, Y,h) = j e'^''^''^''^^^^^-^^xv''^\o{v)dT'dadr] . 

By fl2.76p and fl2.8ip . Proposition 12.141 is clearly a consequence of the following estimate: 

(2.83) Yl -^\lN,2{T,X,Y,h)\ < C~X-'/' 



4Ar/T6 [1/4,2] " 

with C independent of X G [0, 1], T g]0, a"^/^], Y e R, a e [/i", oq] and A G [Aq, oo[ with 
ao small and Aq large. 

Lemma 2.23. For all k there exist Ck independent of X e [0,1], T G]0,a"^/^], y G M, 
a G [/i", ao] and A G [Aq, oo[ such that 

(2.84) J] -^\^^4T,X,Y,h)\<aX~'. 

4Ar/Te [1/4,2] 

Proo/. Recall that T' = -2^^^ + a + a/i^. Let us define the functions (see fl2.28p ) 

(2.85) y = - <y)Hi + \{a' - f?) + 4iV(l - -B'{i^X{l + /i^)t))i/2 



3 



r = 2v/p + a/i2(a - /i + 2iVyrT7^(l - ^5'(r/A(l + /i^)^)) . 
There exists a universal constant Co, with $ = ?7(F + 4>aNri\)^ such that 

(2.86) \x -x\ + \Y -y\ + \T-r\< Co(|$;i + + \K\ + ■ 

The left hand side of f l2.86p does not depend of z. Since G = ri{Y + (t>a^N,r]\) is a critical 
value of $ with respect to z, we get 

(2.87) |X - A"! + |F - 3^1 + |T - n < Co(|G;| + + |G;|) . 

Consider any given (X, Y, T). Then for X ^ A/i(X, F, T), we have |x— + 3^| + |t— T| 7^ 
for all values of fi, a, 77, A, all |x — X| < 1, all — F | < 1, all \t — T\ < 1. This implies 
\X-X\ + \Y-y\ + \T-r\ > 1. Therefore, for all values of /x, a, r/. A, and X ^ M (X, F, T) 
we get 
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From fl2.79p . using 77^^ = \d\, N/X^ G 0{h^) and the fact that any function of type iVOs 
is in 0{N~'^) , we get that all the derivatives of G with respect to rj^T'^a are uniformly 
bounded. By integration by part in (12.82^ we thus get |7Ar.2| < Ck\~^ for all k, with Ck 
independent of X G [0,1], T e]Q,a-^% F G M, and ^ M(X,r,T) . The proof of 
Lemma [2.231 is complete. □ 

From Lemma [2]23l and since |A/i(X, F, T)| is uniformly bounded by lemma 12.1 8[ we get 
that fl2.83p will be a consequence of 

(2.88) ViV with m/T G [1/4,2] , I e'^^'^-^-^xdT'da\ < C\~^/\ 

yN J 

Here and after, we denote by C any constant which is independent of > 1, X G [0, 1], 
T g]0, a"""^/^], a G [/i°, ao] and A G [Aq, C)o[ with small and Aq large. 
Observe that we can now replace the phase '^a,N,\ by the phase 

(2.89) i;a,N,x = {X- Fo)a + a^/S + HoT' + MT') 

+ 2iV(l+U)V. (-(^ + «^o)^/^ + T'/2r - ^(T72 + a)3 + aNO. 

since by (12.79^ the difference '^a,N,\ — iipa,N,x + 5'i/A^) does not depend on T', a, and e^^'^^^ 
is a classical symbol of order in A. We set x = e*^^/^x and recall h = a?/"^ /\ . Then 
x(T', cr; T, a, N] A) is a classical symbol of degree in A > Aq, compactly supported in 
(T', a) close to {0} x [—1, 1]. Moreover, for all a, there exists Ca independent of a, N and 
T G [1/4,2], such that sup(r, ,^;^) l^{T',a)X| < C'a- 

Lemma 2.24. There exists C such that for all N > \^/^ , 

(2.90) -^1 j e'^^'^'^'^xdT'dal < C\-^^\ 



Proof. It is sufficient to prove that, for all N > \^^^, 

(2.91) I j e'^^--«'^xdT'da\ < CA'^/^ . 

Set X - Fq = -AA"2/3^ifo = -5A^2/3 perform a change of variable in ( ICT]) : T' 
\~^/'^x, a = X'^^^y . We are reduced to proving 



(2.92) I j e'^x{^-^'^x,\-^''^y, ...)dxdy\ < C 
with a phase G of the following form 

(2.93) G = -Ay + y^^ - Bx + XU>^-'/'x) + ^j^^^^f^p^y/^ ivi^ + <^F^f" + ^^f 
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Then (12.921) is an oscillatory integral over a domain of integration of size A^^^. Parameters 
-Fq, Go, X^^^/N are bounded, and the main point is to prove that the constant C is uniform in 
{A,B) = {r cose, r sine) with r < cqA^/^. Recall X^PaiX'^^^x) = ^^^^^^^{l+0{ax^X-~^/^)). 
One has 

d^G = -B + X'/'i>',iX-'/'x) + Oiix, y)) + 0(1) 

dyG = -A + y' + 0{{x,y)) + 0{l). 

Moreover, since x is compactly supported in (T', a) one has, with Ga independent of 
T, a, N, X, 

sup ^,)x(A-^/'x, A-i/^y, ...)| < C,(l + \x\ + . 

Therefore, for any Tq, the oscillatory integral (I2.92p is clearly bounded for < r < Tq 
(integrate by part for large {x,y)). 

For r G [ro,coA^/^], we rescale variables {x,y) = r^^'^{x' ,y'), and we set G = r^l'^G' and 
x'(a;', y', .••) = r"^/^A~^/^l/', .••)• Observe that since r^l'^X~^l'^ is bounded, we 

still have decay estimates 



sup <C,(l + k'| + |y|)■ 

(x^y) 



-|a| 



We have to prove 

(2.95) r| y e^^'^'^'x't^x'rfi/'l < C. 

First, the critical points of G' satisfy 

d,,G' = -cose + r-iA2/3^;^(rV2A-V3a.') + r-^/^0{{x', y')) + r-^0{l) 

^^■^^^ dyG = -sme + y'^ + r-i/2(3((x', y')) + r-^O(l) . 

Let F{u) denote a smooth function near n = 0, we then have 

r-^X^'^'Sr^'^X~^'^x') = --^^—-{l + arX-^'\'^F{a^'\^/-'X~^/\')) . 

8(1 + aY'^ 

By (I2.96p . the contribution of large {x',y') to (I2.95P is 0(r~°°), and we may localize the 
integral on a compact set in (x', y'). For | sin0| > 0.1, and tq large, we have two distinct non 
degenerate critical points y'^ = ±| sin^|^/^ + 0(r~^/^) in y' with critical values G'_^{x', ...) 
and thus we get by stationary phase 

(2.07, ./e--V<^W . ^"(je-"o'.^^ . 

Moreover, one has 

9,,G;(x',...) = a,.,G"(x',i/;,...) = - cose + r-^X^'^'^{r^'^X-^'^x') + 0{r-^l^) 
and this imply > Cq > 0. Thus, we get 
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Therefore, ()2.95p holds true, and in fact, we have the following better estimate in the range 
r > 1: 

(2.98) r| / e"''''^'^dx'dy'\ < Cr~^'^ . 



If sin 6' is close to 0, then we first perform the stationary phase in x', and we use the same 
arguments. This completes the proof of Lemma 12.241 □ 

Therefore, we can now assume that = A > 1, and take A as our new large pa- 
rameter. We set X — Fq = —pN'"^, Ho = and we change variables in fl2.88p : 
T' = -2x/N, a = -y/N. We will prove 



(2.99) I j e'^^^x{x/N,y/N,...)dxdy\ < CA'^/^ 

with a phase Qa which takes the following form 



(2.100) ga=py- yVS + qx + N^a{-2x/N) + ^^^-^^(y(l + aF^f/^ + xf 

Observe that, for large N, (12.991) implies a better estimate that f l2.83p : more precisely, 
(I2.99P is equivalent to 

(2.101) -^1 je'^'^-'«'^xdT'da\<CN~^/^\-^/^. 

The above estimate is of course compatible with f l2.90p for ^ A^^/^. 

Recall, see (ET]), that 'ipaiT') = ^^^^^{l + 0{aT'^)). From fl2:93|) we get 



d.Ga =q-x' + 7T-%TT7^(y(l + aFo)'/' + x) + -^,{x + yf + aO{{x, yf) 
(2.102) (l + «i^o)^/2 4Ar2 



dyGa =P~y^ + Go{y{l + aFoY'^ + x) + ^(x + y)2 + aO((x, yf) 



and 



•^L^a = - 2X 



(2-103) dl^G^ =Go + ^,ix + y) + aO{{x, y)) 

dlyGa = -2y + Go{l + aFo)i/2 + + ^) + ^o{{x, y)) . 
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Thus we get the value of the hessian Hn^x, y; T, a) 

(2.104) HNix, y; f, a) = det (^|'^|^ 

= -2Go(x + y)+4xy-^{x + yf + aO((x, y)) . 
Lemma 2.25. There exists Tq and C such that for all {p,q) with \ {p,q)\ > Tq 

(2.105) I / e'^^-'xix/N, y/N, ...)dxdy\ < CA"^/^ . 



Proof. Set {q,p) = (r cos 6*, r sin 6*) with r > Tq large. Let x ^ C*^(|(x,?/)| < c) with small 
c and X = 1 near 0. Then by (12.1021) we get for all k by integration by part in (x, y) 

^'^^^X{r~^^\x,y))x{x/N,y/N, ...)dxdy\ < Cfcr'^A"*^ . 

For the remaining term, we perform the change of variable {x,y) = r^/^(x',y') and we set 
Q'^ = r~^^'^Qa- It remains to prove 

(2.106) |r [ e'^^'^'^^'^il - x)ix' ,y')xir'^^x' /N,r^/%' /N, ...)dx'dy'\ < CA-^^^ . 



Observe that since (1 — x){x'^ u') = near 0, (1 — x){x'^ v') = 1 for \{x'^ y')\ ^ c, and since 
x{u,v, ...) is compactly supported in {u,v), we still have 

s«P(x'y)l^r-'y)(l - xKx\y')x{r'/'x'/N,r'/'y'/N, ...)| < C„(l + \x'\ + |y'|)-l"l . 
The phase Q'^ is of the form 

= sin Oy'-^ + cos 9x'-^ + +-^{x' + y^ + r-'/'0{{x', y'f) + aO{{x', y'f) 

where 0{{x' ,y'Y) means any function of the form x'^y"" f{r^^'^x' /N,r^^^y' /N, a, N) with / 
smooth uniformly in a, N and j + 1 = k. Thus for small a and large tq, we may localize the 
integral (12.1061) to a compact set in {x',y') (integrate by part). The hessian of G'^ is equal 
to Ax'y' — -^{x' + y'Y + 0(r~^/^ + a). Thus for N > 2, a small and tq large, the set on 
which the hessian vanishes defines a smooth curve F outside (x', y') = (0, 0), which is close 
to the union of the two lines c{x' + y') ± (x' — y') = 0, = ^^^^ E [3/4, 1]. Moreover, one 
has 

d^^g', =cose- x'^ + 7^(x' + yT + 0{r-'/' + a) 
(2.107) 4iV 



dy,g'^ =sm9- y'' + ^(x' + y'f + ^(r-V^ 



The contribution of points {x',y') outside F to the left hand side of (I2.106P is estimated 
by 0(r~^/^A~^) by the usual phase stationary theorem. To estimate the contribution of 
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points {x',y') close to F, we use lemma [2.201 For any value of 9, one gets easily that the 
hypothesis of part (a) of lemma [2.201 holds true, and this gives the estimate 

,-^/^AGi(l - x)ix',y'Mr'/\'/N,r'/y/N, ...)dx'dy'\ < Cir^'/^A)-'/' = r'^/^A-^/^ 

which gives the bound Cr~^^^A~^^^ on the right hand side of (12.1061) . The proof of Lemma 
12.25! is complete. □ 

We can now assume \{p, q)\ < tq. There exists c > independent of > 2 such that 

(2.108) V(a;, y) G R\ \x^ _ _l_(:r + ^)2| + |^2 _ ^ ^)2| > ^^^2 ^ _ 

Thus by integration by part, (12.1021) yields that large values of {x,y) gives a contribution 
0(A~°°) to the integral (I2.99P . We can then replace xi^/^, u/N, ...) by a symbol xi.^^ •••) 
compactly supported in the ball B = \ {x,y)\ < R with R large. We are left to prove 

(2.109) ^ |/e-,<^.<^.l<CA-"^ 

Uniformly in A^ > 1 and T near [1/2,3/2] and {x,y) near B, we have 
Ga =Go + 0(a) 

(2.110) f 1 

Qo =qx - xV3 + py- y'/S + -{x + yf + Y^i"" + ■ 

Note that the hessian of Qa is Tia = — 2T(x + y) + ^xy — ^{x + y)"^ + 0(a). Therefore the 
set Za = {{x,y)]'Ha{x,y) = 0} is, for small a, a smooth curve in B which is close to the 
parabola —2T{x + y) — {x — yY = for A^ = 1, and close to the hyperbole —2T{x + y) + 
4:xy — ^{x + yY = for N > 2. Moreover, 

dxGa =q - y, a); X{x, y, a) = x^ - f{x + y) - 7^(x + yf + 0{a) 

(2.111) Y 

dyGa =P - yix, y, a); y{x, y, a) = y'^ - f{x + y) - -^{^ + vf + 0(a) 

It remains to use Lemma [2.201 near any point {q,p), \{p,q)\ < ''"o- If (q^p) is not in the 
image of Za by the map {X,y), then near {q,p), the estimate (12.1091) holds true with a 
factor CA~^ on the right hand side by the usual stationary phase theorem. If {q,p) is in 
the image of Za by the map {X, y), but {q,p) ^ (0, 0), then one easily verifies that part (a) 
of lemma 12.201 applies, and this gives near {q,p) the estimate (I2.109P with a factor CA^^/^ 
on the right hand side. Finally, near {q,p) = (0,0), one has {x,y) near (0,0), and one 
easily verifies that part (b) of lemma [2.201 applies, and therefore (I2.109P holds true. 

This concludes the proof of Proposition I2.14[ 

Remark 2.26. Since the symbol x of degree is elliptic at {x,y) = (0,0), the estimate 
(I2.109P is optimal. To see this point, it is sufficient to apply part (b) of lemma 12.20! at 
(p,g) = (0,0). Observe that by f[2T8D . (p, g) = ( 0, 0) is eq uivalent to X = Fq = 1 and 
T = + a, i.e equivalent to x = a, t = 4:N^y a(l + a) which are precisely the times 
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where a swallowtail occurs in the wave front set of the Green function. This, and f l2.10ip . 
proves Theorem II. 4[ 

2.6.5. Proof of Proposition \2AR In order to prove Proposition 12.151 we use as before the 
splitting 7i = 7i^i + 71^2- First, we prove 

(2.112) |7i,i(T,X,r,/.)| <C(27r/.)-2((-^)V2 + /,V3). 

Going through the proof of Proposition I2.13[ we notice only one difference between the 
case = 1 and N > 2: namely for T g]0, Tq], when we estimate J e*'^*^'+'+G+^+(2;; aX)dz, 
we may have a critical point associated to a very large value of z. Let x{^) ^ ^'^(j^i, 00 [) 
with Zi large, and set 



(2.113) J = J e'^^''+'+Q+^+{z;aX)x{z)dz. 
We will prove 

which clearly gives the first term on the right of the inequality fl2.112p . One has 

5.^1,+,+ =^(1 + az)-'/' - ^(1 +X) + 0{z"^'') 

(2.114) ^ ^ 

=^(1 + az)-^" + ^(1 + X) + 0{z-''') . 

Therefore, to get a large critical point Zc-, T must be small (recall T < Tq means t < 
a^/^To). One has then Zc^^"^ ~ T and from fl2.114p we get a|$i,+,+ (^J ^ T^- Recall that 
9_|__+(z, a. A) is a classical symbol in z of degree —1/2, thus T~^s^/^9+^+(T~^s, a. A) is a 
symbol of degree in s > sq > 0, uniformly in T g]0,To]. Therefore, if we perform the 
change of variable z = T~'^s in (12.1130 . we get \ J\ < CA~^/^T~^/^ by stationary phase. 

It remains to prove that the inequality (I2.26P holds true for 71 2. The only place where 
N >2 gets used in the proof of Proposition 12 . 141 is Lemma [2l25] and inequality (12.1080 . But 
for = 1, since x(a;,|/, ...) is compactly supported in {x,y), we do not need the inequality 
f l2.108p . Moreover, we get from (12.1020 that the phase Qa has no critical points on the 
support of X ioi |(p, q)\ > tq if tq is large, and this implies for |(p, q)\ > vq 

(2.115) I j e'^^''x{x,y,--)dxdy\ e 0(A"°°) . 

In fact. Lemma [2.251 is telling us that the constant C in the right of (I2.105P is uniform for 
large A^. The proof of Proposition 12.151 is now complete. 

3. Parametrix for < a < h^^"^ 

We will write the initial data with the help of gallery modes, which we first describe in 
connection with the spectral analysis of our Laplace operator. We describe the correspond- 
ing solutions of the wave operator. We then estimate their L°°(f2) norm for tangent initial 
directions by using Sobolev embedding, taking advantage of the size of the Fourier support. 
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We deal with the non-tangent initial directions by constructing a crude parametrix, relying 
partly on gallery modes and the asymptotics of the Airy function. 

3.1. Whispering gallery modes. Let = {{x,y) E M?\x > 0,?/ G M} denote the half- 
space with the Laplacian given hj Ad = + {1 + x)dy with Dirichlet boundary 
condition on dQ. Taking the Fourier transform in the y-variable gives 

(3.1) - A^,, = -dl + (1 + x)t]^. 

For ?7 7^ 0, —Ad,ti is a self-adjoint, positive operator on L^(M+) with compact resolvent. 
Indeed, the potential V{x,ri) = (1 -|- x)?7^ is bounded from below, it is continuous and 
lima;_>.oo V{x,ri) = oo. Thus one can consider the form associated to —d^ + V{x,ri), 

Q(u) = I \dM'^+V{x,7])\v\^dx, D{Q) = H^(R+)n{v e l\r+), (1+x)^/\ g i:^(M+))}, 

Jx>0 

which is clearly symmetric, closed and bounded from below by a positive constant c. If 
c ^ 1 is chosen such that —Ao,r] + c is invertible, then {—Ao,ri + c)~^ sends L^(]R_|_) in 
D{Q) and we deduce that (— A/) + c)"^ is also a (self-adjoint) compact operator. The 
last assertion follows from the compact inclusion 

D{Q) = {v\d^v, (1 + xy/'^v G L2(R+), v{0) = 0}^ L\R+). 

We deduce that there exists a base of eigenfunctions Vk of — A^^^ associated to a sequence of 
eigenvalues Xk{r]) ~^ cxo. From — A^,,!; = Xv we obtain d^v = (77^ — A -|- xrf)v, w(0,?7) = 
and after a suitable change of variables we find that a basis of L^([0,cxd[) orthonormal 
eigenfunctions is given by eigenfunctions 

(3.2) ek{x,r]) = fk-g^Aiirfix - Uk), 

where {—u:k)k denote the zeros of Airy's function in decreasing order and where fk are 
constants so that ||efc(., ?7)||l2([o,oo[) = 1 for every /c > 1, and the fk remain in a fixed 
compact subset of ]0, oo[. The corresponding eigenvalues are Xkijf) = tf + (^kV^ ■ 

Remark 3.1. Let 6x=a denote the Dirac distribution on M_|_, a > 0, then it reads as follows: 

(3.3) 5a;=a = ^ efc(x, r])ek{a, 77). 

fc>i 

We define the gallery modes as follows: 
Definition 3.2. For a; > let Ek{^l) be the closure in L'^{Q) of 

(3.4) 1^ j e'y'^ek{x,v)Mv)dv,^k e 5(M)}, 
where S{M.) is the Schwartz space of rapidly decreasing functions, 

5(M) = {/ G C°°(M), < 00 Va,/3 G n}. 

For k fixed, a function in Ek{^l) is called whispering gallery mode. 
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We have the following result (see |6]): 

Theorem 3.3. We have the orthogonal decomposition {L^{Q), Ad) = @j_Ek{Q), where 
Ek{^l) denotes the space of gallery modes associated to the kth zero of the Airy function 
Ai and where Ad = + {1 + x)dy with Dirichlet boundary condition on dQ. 

Proof. Indeed, from [6] [Section 2.2] one can easily see that {Ek{fl))k are closed, orthogonal 
and that UkE^i^l) is a total family (i.e. that the vector space spanned by UkE/^^Q) is dense 
in L\n)). D 

Let tpj G C^(]0,oo[) as in remark [1.101 Using remark IXT| for h G (0, 1], we write the 
initial data, localized at frequency j^, as follows 

(3.5) iJ2{hy^ -Ag)'ipi{hDy)6^=a,y=o = 

^^rh / ^^^''^2(?7-\/l + uJk{h/r]y/'-^)tJji{r])ek{a,r]/h)ek{x,r]/h)dr] . 

l<k 

Observe that in the sum over k, by remark [1.101 we may assume k < eh~^ with e small. 
From ( 13. 5 p we get 

(3.6) Ua,hit,x,y) = e-'''/^^Mhy^^g)MhDy)6^=a,y=o 
i<k ^ ^ 

Our goal is to prove the following proposition. 

Proposition 3.4. There exists C such that for every h g]0,1], every < a < h^^"^ and 
every t E [—1, 1], the following holds true 

(3.7) ||l,.<aMa,;,(t,x,y)|Uo. <C/i-2min(l,/ii/^ + (|^)i/=^). 

This proposition will be proved in the next two sections. Proposition 13.41 clearly implies 
theorem 11.31 for a < h^/"^ . By time symmetry, we may restrict ourselves to positive times 
t G [0,1]. Notice that the proof for the wave propagator ey:.\;){+^t^J —Ag) is exactly the 
same as the sign plays no role whatsoever. 

3.2. Tangential initial directions. In this section, we make use of the Sobolev embed- 
ding properties related to the orthogonal basis (e/o). 

Lemma 3.5. There exists Co such that for L > 1 the following holds true 

(3.8) supf V k-^^^At'ib-iOk)) <CoL^^^. 



l<fc<L 

Proof. From |Az(x)| < C(l + |a;|)~'^/^, we get 



l<k<L l<k<L ' 
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From cjfc — k"^^^, we get easily with C independent of L and D large enough 

sup Jib) < CL^/\ sup J{b) < CL^/^ . 

Thus we may assume b = with b' e [0,D]. Since cok = k'^^^g{k) with g being an 

elliptic symbol of degree 0, we are left to prove that 

H^) = L-"^ E (^/^)"'^'l + LV3|a:-(A:/L)2/3|i/2 

l<k<L ' \ I I \ 



satisfies sup2.g]u /(x) < CqL^^^. Since we can split [0, 1] into a finite union of intervals on 
which the function ^_^i^i_/l\^^^^2/i\i/2 is monotone, and since each term in the sum is bounded 
by 1, we get 

/(.) < Cte + L^' [ —£^—^^it < Cte + [ ^^-^d. . 

and the proof of Lemma 13.51 is complete. □ 

Let Ua,h,<L be the function defined by (13. 6p with the sum restricted to k < L. From 
(13. 2p . lemma IX^ and Cauchy-Schwarz inequality, one gets 

(3.9) \K,h,<L{t,x,y)\\L^ < C,h-'h'/^L'/\ 

Taking L = C/h, we get that proposition 13.41 holds true for |t| < h. With L = h~^^^ 
OT L = l/\t\, one sees also that (13.71) holds true for Ua^h,<L- Thus we are reduce to 
prove that (13.71) holds true for Ua^h,>L which is defined by the sum over k > L with 
L > D raaxi^h"^^^, Vl^l) "with a large constant D . 

3.3. Non-tangential initial directions. In this section, we denote by Va^h{t,x,y) the 
function defined for /i < t < 1 by (13. 6p with the sum restricted to L<k<e/ h, L> 
D max{h~^^^, 1/t), D > large and e > small. For each value of k, we set 

(3.10) X = tujkh-'/', /i = ^^. 

From Uk — k'^^^, k > 1/t, and t > h, one has for some c > A > c; thus we will take A as 
our large parameter. However, the parameter /i just satisfies 

(3.11) < ^ < ^ mm{t^/'\ h^'^^) 

and thus may be small or arbitrary large. Observe that for D large enough, for k > Dh~^/^ 
and < X < a < h^^"^ one has 

for all rj in the support of ipi. Therefore we can use the asymptotic expansion of the Airy 
function from section I2l2l with u = e™^^. 
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M{0 = J]a;±^e^i'(-^)^^^-C)-^/^^±(-C) 
± 

which is vahd for — ^ > 1, with —( = uJk — 'rf/^h~'^/'^x > co'fc/2 > ^ > 1 since uji ~ 2, 33. 
We thus get from (13.61) 

(3.12) va^,= ^'^• = ^E /^^'"^'"^f'^^^ 

L<k<e/h ±,± 

The phases ^^'^ and the symbols a^'^ of Wk read as follows, with the notation z = 

h^/^UkT]-^/^ > 2a: 

(3.13) hX^p^it, X, y, T],a) = ri{y- tVTT^ ^)'^' 



(3.14) a^Jx,r],a,h) = h-'/'r]Mv)MvVTT^)-^cu^uj^ 



p2 
_fc_ 

'^1/3^ 

X (^ - xr'/\z - a)^^/^^±(r/2/3/i-2/3(^ _ a;))^^(^2/3/,-2/3^^ _ 

One has Srjdrf = —2zdz and for < x < a < 2z, 

|(za,)^az - xy"')\ < C,z-'" < C'^ihk)-'/' ; 
moreover, \E'-|- are classical symbols of degree at infinity and 

\r^y^h-'/%z-x)\>Uk/2>Ch-'/\ 

since k > L > Dh^^^^. Therefore we get from f l3.14p that for all j, there exists Cj 
independent of h, k, a, x, r] such that 

(3.15) \d'^ap^ix,r],a,h)\ < C,{hk)-^/\ 

Proposition 3.6. For e small, there exists C independent of a E {0,h^^'^], t E [h,l], 
X E [0,a], y G M and k E [L,e/h] such that the following holds true 

(3.16) I J e^"*^'Vj'^dr/| < Cikky^/'' X^'/^ . 
Observe that from f l3.12p and the definition (13.101) of A, ( I3.16P implies 

\\t.<aVaAt,^,y)\\Loo < ch~' Yl {hky^iH~^i^h^'''k~^i'' 

k<e/h 

(3.17) , , 

= C/i-2(-)V3/,i/9( Y k-'/') < C'h-\-f/^ 

k<e/h 

and therefore proposition 13.41 holds true for Va,h- 
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Proof. Since from f l3.15p the {hkY^^a^'^ are classical symbols of degree compactly sup- 
ported in rj, we intend to apply the stationary phase to an integral of the form 



±,± 



(3.18) J = J e'^'^k' gdr] 

with g a classical symbol of degree compactly supported in rj. We have to prove uniformly 
with respect to the parameters the inequality 

(3.19) |J|<CA-^/^ 
Derivative of the phase with respect to rj yields 

hXd,<Pp^ = y - ti±ii ± lx{z - xfl' ± \a{z - af'\ 

yjl + Z O 3 

where the two ± signs are independent from each other (thus, we have 4 cases to consider). 
Let (5 = ^ e [0, 1], a = ^^2^1^ and s = r/"^/^ g [sq, Si]. Since D is large, one has a G [0, Cq] 

with Co such that r/"^/^ = s > Sq > 2co on the support of i'lir]). Let X = ^J~^2/3 and define 
the function g{z) by 

Then the derivative of the phase is equal to 

(3.20) drj(j)p^ = X- sg{h^'^ujks) + ^/i^^'^, ^^'^ = ±5{s - 6af'^ ±{s- af'^. 
We now study critical points. We take s = r]^"^/^ as variable and we get 

9.9,0^^ = -{g{z) + zg'{z)) + U±5{s- 6a)-'/' ± {s - a)~'/'" 

(3 21) 

dld,4>i^ = -0'u,{2g\z) + zg\z)) " f ( ± ^(s - 5a)-'/' ± 
Lemma 3.7. For e small enough, there exists c > independent of k < e/h such that 

Proof. One has (s — a)-'/' > 6{s — Sa)-'/'; for e small, z = h'/'tOkS is small and thus 
g{z) + zg'{z) is close to |. Thus we get 

\dsdr,cpp-\ > 1/10. 

The derivative dsdr^(t)j^ ' may vanish but in case \dsdrj(l)^'^\ < 1/100, the first line of (^23) 
implies 

^(s_tt)-i/2>o.05. 



[s — a] 
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The second line of (13.211) then gives a positive lower bound on \d'^dri(f)p~^\. It remains to 
study 0^'"''. For any function /, one has 

/•q(1-(5) 

(3.22) f{s -a)- 5f{s - 5a) = (1 - 5)f{s -5a)- / f'{s -5a- t)dt . 

Jo 

Taking f{t) = t~^^'^, we thus find that 

'+1 < 1/100 ^ ^(1 - 5) > c > . 

If one applies (I3.22p with f{t) = t~'^^'^, we then find that for e small, the second line of 
13.211 implies I^^S^^^'"^! > c/2. The proof of Lemma \377\ is complete. □ 

From lemma 13.71 and I2.19[ we get that proposition 13.61 holds true in the case where the 
parameter /i is bounded, since in that case all the derivatives of order > 2 of the phase 
0^'^ are bounded. It remains to study the case where fi is large. 

In cases (+, +) or (— , — ), and fi large, we can take as large parameter A = A/i. Since 
(s — a)"^/^ + 5{s — 5a)~^^'^ > c > 0, we get in that case that (I3.16P holds true with a better 
factor (/i/c)~^/^A~^/^ on the right hand side. 

It remains to study the cases (+, — ) and (— , +) for fi large. But in that cases, we can 
use (13.221) : therefore, if /x(l — 5) is bounded, all the derivatives of order > 2 of the phase 
(pp^ are bounded, and therefore from Lemmata 13.71 and 12.191 we get that proposition 13.61 
holds true. 

Finally, in the cases (+, — ) and (— , +) and /i(l — 5) large, we can take as large parameter 
A' = A/i(l — 5), and since by (13.221) one has 

\{s - a)-i/2 _ _ 5ct)-i/2| > c(i _ S) , 

with c > 0, we get in that case that (13.161) holds true with a better factor (/i/c)~^/'^A'~^/^ 
on the right hand side. 

The proof of Proposition 13.61 is now complete. □ 

This concludes the proof of Proposition 13.41 

4. Dimension d>3 

Let d > 3 and fi^ = {(x, y) x M*^"^} with Laplace operator = + {1 + x) Ay. 

The normal variable is still denoted x > 0, and the boundary is still defined by the condition 
X = 0. The proof of Theorems 11.31 and 11.41 follows exactly the same lines as in the 2d case, 
both in the case a < h^^'^ and a ^ h^^"^. 

4.1. Parametrix for a 3> /i^^^. In higher dimensions the parametrix construction is iden- 
tical to the one in the two dimensional case. We set h = h/\r]\, and we define v(t,x,y, h) 
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with y e M'^"^ by 

v{t,x,y,h)= ^ VNit,x,y,h) 

with the same UN{t,x, h/\ri\) as before. We take polar coordinates in 77 G M''""^, rj = \r]\u. 
We thus get 

(4.2) v^{t,x,y,h) = J ( J e''-^'^y-'^Uu)u^{t,x,h/\7^\M'xo{\r]\)d\vl 

In the above formula, apart from the harmless factor \r]\'^ instead of |?7|, we have a superpo- 
sition with respect to a; G S'"'"^ of functions of the same type as before, which are evaluated 
dX z = y.oj. We shall use the following lemma. 

Lemma 4.1. Let ipj G C^(]0, oo[). There exists cq > such that for every a g]0, 1] and 
every t G [h, 1] 

(4.3) /l'|V^l(/lv/^)^2(/i|/^,|)e^**^<5.=a,j,=o|Loo(.<a,|j;|<cot) £ 0((^)°°) 

Proof. We may and will assume a < 2t; In fact, for t < a/2, by finite speed of propagation, 
the singular support of e^'^^^~'^'^6x=a,y=o has not reached the boundary x = 0, and then (14. 3 p 
is a simple consequence of propagation of singularities in the interior (see the argument 
below). Let T G [h, 1] be given; perform the change of variable t = Ts,x = TX,y = TY, 
and set /r(s,X, F) = f(Ts,TX,TY). Then one has 

{AJ)t = T^^PrfT, Pt = 4 + (1 + TX)Ay 
Set h = h/T < 1. One has for any ip the identity {;ilj{hDt^x,y)f)T = V^(^-Ds,x,y)/r, and 
therefore (14.31) is equivalent to the estimate at time s = 1 



(4.4) |^i(/iVi¥)V^2(/i|I5y|)e^*^5x=a/T,,=o|L-(x<a/T,|y|<co) e 0(n°^) . 

Observe that b = a/T < 2 is bounded. Since '?/'2(^|-Dy |) commutes with the flow e^*^^~^, 
using the Melrose- Sjostrand theorem on propagation of singularities at the boundary |10] . 
we just need to verify the following : There exists Cq > such that for any T G [0, 1] 
and any optical ray s — > p{s) associated to the symbol + (1 + TX)ri'^ starting at t = 
from p(0) = {X = h,Y = 0;^o,^o) with + (1 + Th)ril = 1 and |r/o| > Ci > 0, one 
has |y(p(l))| > 4co. But on the generalized bicharacteristic flow, one has dsTj = and 
dgY = 2?7(1 + TX{s)) and therefore Y{s) = riQ{g{s)) with g{s) > 2s and the result is 
obvious. Observe that the cutoff by ■ip2{,h\Dy\) is essential to get the lower bound on |?7o|. 
The proof of Lemma [4.11 is complete. □ 

In order to prove our dispersive estimates, we may assume h < t < 1, and therefore by 
lemma WTLi we may also assume \y\ > cot > cqH. Classical stationary phase in G S'^~^ 
gives 
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/ A ^\ f ,-Mi'„,,,'i , / h N , '\v\\y\ , h , -i\v\ly\ , h 

(4.5) e' i^^y-'^'du = {----) ^ {e a +{----)+ e ^ 

J my\ my\ \v\\y\ 

where a± are classical symbols of degree in the small parameter j;^- Inserting f l4.5p in 

(14. 2p . and since for \y\ > c^t and \'q\ G [|,8] one has (j^^)^ ^ C{h/t)^, we easily see 
that the proof of theorem II .31 and 11.41 follows exactly like in the 2d case. 

4.2. Case a < h^^'^. Indeed, the dispersive estimates follow once we notice that definition 
13.21 and theorem 13.31 extend to the d dimensional domain fi^. It is enough to define for 
X > 0, Ekiytd) to be the closure in L'^iVL^) of 

(4.6) {-(2^^ j e^^'''^^^(l^l^^-^fc)^(^)^^'^e5(R'^-i)}, 
where 5(M'^~^) is the Schwartz space of rapidly decreasing functions, 

SiR'-^) = [fe C°°(M^-i)|||^"Z}^/|Uo.(M.-i) < oo Va,/3 G n'^^]. 

Theorem 4.2. We have the orthogonal decomposition {L^{Qd), ^d) = ®±Ek{^ld), where 
Ek{^ld) denotes the space of gallery modes associated to the kth zero of the Airy function 
Ai and where = (9^ + (1 + x)Ay with Dirichlet boundary condition on dQd- 

Therefore, by lemma HIT] and (14.51) . the proof of our main theorems follows exactly like 
in the 2d case. 

Appendix A. The energy critical nonlinear wave equation 
We consider the equation 

(A.l) agu+\u\^u = 

with data (uo,ui) G HQ{^ld) x L^(f2d), with 3 < < 6. When the domain is W^, there 
is a long line of seminal works regarding this model, which may be one of the simplest 
model of a critical wave equation. To our knowledge, the fist work to address the energy 
setting (as opposed to C°°) is [12], where low dimensions are dealt with, using only the 
oldest Strichartz estimates (time and space exponents are equal). Higher dimensions (e.g. 
d > 7) have their own set of difficulties, mostly related to the low power nonlinearity 
(1 + 4/((i — 2) < 2) and the subsequent failure of its derivative with respect to u to be 
Lipschitz. All these technical annoyances may be solved one way or another, but are out 
of the scope of the present paper. 

Hence we contend ourselves with the low dimensions. There are essentially two things 
to be checked: 

• we have a "good" local Cauchy theory, providing energy class solutions ; 
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• this local Cauchy theory may be tweaked as to insert (a small power of) the potential 
energy of the solution in the nonlinear estimates, so that we can then perform the 
non concentration argument from [3] and extend our solutions globally in time. 
Remark that the potential energy is which corresponds to the critical 

nature of the equation, as by Sobolev embedding, M- 

We refer to [7], [8] for details on how to deal with fractional derivatives, Besov spaces on 
domains and product-type estimates (alternatively, one may proceed with interpolation as 
in 0]). 

Remark A.l. Note that in proving Theorem 11.61 from Theorem II. 3^ one needs, for p > 2, 
an embedding 5^'^ C on domains, which may be proved directly or follows from a 
Mikhlin-Hormander multiplier theorem from Alexopoulos (see [7] and references therein). 

Having these tools at hand, we may proceed exactly as in M"', provided we have the right 
set of exponents. 

• Case (i = 3: Theorem II .61 allows for the Strichartz triplet (g = 4, r = 12, /3 = 1) and 
one may proceed like in the case. This was already observed in [2] and allows 
for a streamlined argument when compared to [3]. 

• Case d = 4: Theorem 11.61 allows for the Strichartz triplet (g = 11/5, r = 22/3, /3 = 
1). As by Sobolev embedding we have M- L^, we may write 

< \u\'">\u\'^/'> e L^Ll X L]L'/' C LlLl , 

and we may proceed as in M^. 

• Case d = 5: Theorem 11.61 allows for the Strichartz triplet (g = 2,r = 5,/3 = 1). As 
by Sobolev embedding we have Hq ^ Lx^^^ , we may write 

H'^'u < \u\^\u\'/-' e L]Ll/' X L^Lio C L]lI , 

and we may proceed as in R^. 

• Case d = 6: Theorem 11.61 allows for the Strichartz triplet (g = 2,r = 18/5, /3 = 1). 
By Sobolev embedding -8^05^ L^, we get u G L'^L'^ which provides a local Cauchy 
theory but without the potential energy factor. However we may estimate 

\u\u e LfLl X L'^B'lj^jf^ c , 

which is the dual endpoint Strichartz space. As we may estimate the B.^^^'-^^ norm 
of u in term of Hq and norms, we now have a good local Cauchy theory, suitable 
to globalization in time. 
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